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ABSTRACT 


Over the past several decades, the Horn of Africa (HOA) has experienced 
recurring climate variations, including droughts and floods that have devastated 
the region’s livelihoods and prompted increased investment in strategies to 
minimize the negative effects of climate variations. These preventative strategies 
include the enactment of early warning systems, such as the Famine Early 
Warning System Network, and military commands such as U.S. Africa 
Command. If these organizations are to be successful, they must account for 
the many climate factors that affect Africa, including seasonal climate variations 
and climate change. Thus, skillful long-range forecasts, especially of 
precipitation, have become increasingly valuable in planning the operations of 
these organizations. 

In this study, we focused on assessing the potential for predicting HOA 
precipitation rate (PR) during the October-November rain season at lead times of 
several seasons. We correlated HOA PR and remote climate variables, and 
discovered a strong potential for skillful long-range forecasts of HOA PR using 
sea surface temperatures (SST) near New Zealand, the Philippines, and Namibia 
as predictors. Our forecast methods included deterministic (tercile matching, 
linear regression, optimal climate normals) and probabilistic (composite analysis) 
methods. Our verification metrics showed a definite improvement in forecast skill 
over existing long-range forecasts based on long-term means, and indicated that 
our forecasting methods have the potential to improve the planning of military 
and non-military operations in the HOA. 
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Figure 1. Schematic illustration of anomalies in the Indian and Pacific basin 
Walker circulations, and in low-level atmospheric pressures and 
sea-surface temperatures, during EN and positive lOD phase 
events. These anomalous conditions tend to contribute to above 
normal precipitation in the HOA. These anomalies tend to be most 
pronounced during August-February, and thus tend to have their 
greatest impact on the HOA short rains period in October- 

December. Figure from Hunt (2009).7 

Figure 2. Outline of process used to generate some of the climate forecasts 
used by FEWS NET. Observations (left box) are used to identify 
statistical relationships between climate system variables (middle 
box; e.g., conditions in the Indian Ocean) and food related variables 
(e.g., precipitation in southern Africa). These relationships are 
used to build predictive statistical models that generate forecasts 
for the food security related variables (right box). Figure from Funk 

(2009).9 

Figure 3. Outline of contingency planning and response process used by 
FEWS NET. Climate forecasts are used as inputs to the 
development of potential food insecurity scenarios that are then 
used to develop contingency plans and food security outlooks (e.g.. 

Figure 4). Figure from Funk (2009).9 

Figure 4. Example of a food security outlook product generated by FEWS 
NET. This outlook is for conditions in eastern central Africa and the 
HOA in April-June 2010. Climate forecasts are an essential 
element in the generation of such outlooks. These outlooks are 
used by agricultural, emergency management, and other planners 
to mitigate the food insecurity events and their effects. Figure from 
U.S. Aid, FEWSNET, May 2010; available online at 
http://www.fews.net/Pages/region.aspx?gb=r2&l=en; accessed May 

2010 . 10 

Figure 5. Example of a seasonal precipitation outlook from ICPAC. This 
outlook shows the predicted probability of above normal (AN), near 
normal (NN), and below normal (BN) precipitation in the greater 
Horn of Africa during March-May 2010. The yellow, green, and 
gray shading represents areas expected to experience similar AN, 

NN, and BN precipitation conditions. Yellow: relatively high 
probability of NN or BN precipitation. Green: relatively high 
probability of NN or AN conditions. Gray: equal chances of AN, 
NN, and BN. The vertical boxes with three numbers show the 
probabilities for each category, with AN at the top of the box, NN in 
the middle, and BN at the bottom. Note that in each colored region, 
the three probabilities sum to 100%. Example: in zone III the 
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probabilities are 35% AN; 40% NN; and 25% BN. Figure from 
ICPAC; available online at www.icpac.net/Forecasts; accessed 

March 2010.12 

Figure 6. Examples of seasonal precipitation forecasts from the Kenya 
Meteorological Department for March-May 2010. March-May in 
central east Africa is normally a period of relatively high 
precipitation that is often referred to as the long rains season. Top 
panel: forecast of areas of above and below normal precipitation. 
Bottom left panel: forecast of onset dates for long rains season. 
Bottom right panel: forecast of cessation dates for long rains 
season. Figure from Kenya Meteorological Department; available 

online at www.meteo.qo.ke/ ; accessed March 2010.14 

Figure 7. Examples of U.S. Africa Command (AFRICOM) activities that 
contribute to U.S. government efforts in Africa. Accurate long- 
range climate forecasts have the potential to improve planning by 
AFRICOM and its African government partners and non¬ 
governmental organization (NGO) partners. Figure from U.S. Africa 

Command (2010).16 

Figure 8. The African region (outlined in yellow) for which the 14th Weather 
Squadron (14 WS) produces its LRFs for the FIOA. The 14 WS 
produces for this region narrative-only long-range forecasts (lead 
times of one to six months) of surface air temperature, precipitation, 
and cloud ceilings. These LRFs are based on composite analysis 
methods employed by the 14 WS, and on LRFs produced by 
civilian operational climate centers (e.g., CPC, IRI, UKMO). The 
locations marked with red dots are ones from which in situ 
observations are collected and for which terminal aerodrome 
forecasts (TAFs) are issued. Figure from personal communication 

with 14 WS.18 

Figure 9. The broad Horn of Africa (HOA) region. The continental scale view 
(left panel) shows the HOA region (red outlined box) in relation to 
all of Africa. The HOA focused view (right panel) shows the 
general HOA region for our study. The darker brown shading 
indicates the countries that are often considered part of the broad 
HOA region. Figure from World Health Organization (WHO); 
available online at: www.who.int/.../hoafrica/en/index.html ; 

accessed May 2010.22 

Figure 10. Long-term mean (LTM) composite of surface precipitation rate (PR; 

mm/day) during October-November (Oct-Nov) for 1970-2009. 
Oct-Nov tends to be a period of relatively high precipitation for the 
general HOA region (see Figures 9 and 11). HOA precipitation in 
this period is: (1) lower than in the western and central tropical 
African and tropical Indian Ocean regions; and (2) higher in these 
months that in the tropical and subtropical regions to the north and 
south of the HOA. Image created using ESRL (2010).23 
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Figure 11. 


Figure 12. 


Figure 13. 


Figure 14. 


Figure 15. 
Figure 16. 


Figure 17. 


Monthly precipitation rate (PR; mm/day; black bars) and standard 
deviation of the precipitation rate (SD; mm/day; white bars) in the 
east African region (5°S-5°N, 35°-46°E) from in situ observations. 
Two rainy seasons occur in east Africa, including the HOA as 
defined in this study: (1) the March to May long rains period; and 
(2) the October to December short rains period. The PR during 
Oct-Nov tends to be relatively high and has a relatively high SD. 
Thus, the Oct-Nov PR is an especially important, but difficult, 
variable to be able to predict at long lead times. The Oct-Nov PR 
in the HOA was the primary target, or predictand, for the LRFs we 
developed and tested in this study. Figure from Behera et al. 2005.. 24 
OAF method equations for forecasting the probability of AN, NN, 
and BN PR (the predictand categories; the terms on the left hand 
sides of the equations) given: (a); the conditional probabilities of the 
PR predictand categories given the occurrence of the specified 
SST predictor categories (the first terms on the right hand sides of 
the equations); and (b) the probability of SST predictor categories 
occurring (the second terms on the right hand sides). See main 
text for more information. The forecast lead time is not specified in 

these equations. Adapted from MetEd (2010).39 

Equations used in our modified CAF process to calculate long- 
range forecasts of PR based on observed AN SST predictor 
conditions at the forecast lead time. The lead time is unspecified in 
these equations. Compare to the standard equations in Figure 12... 41 
Time series of Oct-Nov HOA PR predictand (top panel) and New 
Zealand SST predictor (bottom panel) for 1970-2009. of the bold 
black and red horizontal lines denote the boundaries between the 
AN, NN, and BN categories for the predictand and predictor, 
respectively. For example, PR values are categorized as AN when 
> 3.5 mm/day; NN when < 3.5 mm/day but > 3.1 mm/day, and BN 


when < than 3.1 mm/day.46 

Flow chart showing main datasets and methods used to conduct 

climate analyses and long-range forecasts in this study.48 

The Oct-Nov HOA PR predictand region for our study (black box) 
overlaid on the LTM Oct-Nov PR (color shading; mm/day; cf. 
Figure 10). PR image created using ESRL (2010).49 


Time series of HOA PR (predictand, purple line, mm/day), and El 
Nino-La Nina (ENLN) and Indian Ocean Dipole (lOD) indices (MEI, 
red line, and DM I, green line, respectively) for Oct-Nov 1970-2009. 
Note the similarities in the interannual variability of the predictand 
and the two indices. Indications of teleconnections between HOA 
PR, ENLN, and lOD are highlighted during: (a) a strong LN and 
negative lOD period in Oct-Nov 1974-75; and (b) a strong EN and 
positive lOD period in Oct-Nov 1997. Relationships such as these 
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Figure 18. 


Figure 19. 


Figure 20. 


Figure 21. 


Figure 22. 


suggest the usefulness of evaluating ENLN and lOD indicators as 

potential predictors for HOA PR.51 

Conditional composite anomalies of Oct-Nov precipitation rate (PR; 
mm/day) for the five most extreme AN (wet) and BN (dry) events of 
FlOA PR during Oct-Nov for 1970-2009. The left (right) panels 
show regional (global) scale views. See Chapter II, Section B.2 for 
the five most extreme AN and BN PR event years. The red box 
outlines our HOA predictand region. Note the generally opposite 
patterns in the wet and dry composites, and the correspondence to 
PR patterns associated with ENLN and the lOD. PR anomaly 
patterns such as these suggest possible predictors and linkages 
associated with the climate variations in the HOA. Images created 

using ESRL (2010).52 

Conditional composite anomalies of sea surface temperature (°C) 
for five most extreme wet (left panel) and dry (right panel) HOA PR 
events during Oct-Nov 1970-2009. Note in the wet (dry) 
composite the positive (negative) SSTAs in the WIO and central- 
eastern tropical Pacific and the negative (positive) SSTAs in the 
tropical ElO and western tropical Pacific, consistent with positive 
(negative) lOD and EN (LN) conditions. Images created using 

ESRL (2010).53 

Conditional composite anomalies of 200 hPa GPH (gpm) for five 
most extreme wet (left panel) and dry (right panel) HOA PR events 
during Oct-Nov 1970-2009. Note in the wet (dry) composite the 
positive (negative) GPHAs over Africa and the central-eastern 
tropical Pacific (central and southern Africa and central-eastern 
tropical Pacific), consistent with positive (negative) lOD and EN 
(LN) conditions. These features are indicative of Rossby-Kelvin 
wave responses to anomalous convection over the WIO, MC and 
central-eastern tropical Pacific. Images created using ESRL 

(2010).54 

Conditional composite anomalies of 850 hPa zonal wind (m/s) for 
the five most extreme wet (left panel) and dry (right panel) HOA PR 
events during Oct-Nov 1970-2009. The red arrows represent the 
sense of the wind flow anomalies. Note in the wet (dry) composite 
the anomalous low-level convergence (divergence) over the HOA 
and divergence (convergence) over the MC, consistent with above 
normal PR in the HOA and with positive (negative) lOD and EN 

(LN) conditions. Images created using ESRL (2010).56 

Conditional composite anomalies of OLR (W/m^) for five most 
extreme wet (left panel) and dry (right panel) HOA PR events 
during Oct-Nov 1970-2009. OLR is a good proxy for convection 
and precipitation in the tropics, with positive (negative) OLRAs 
indicating negative (positive) convection anomalies. Note in the 
wet (dry) composite the negative (positive) OLRAs in the HOA, 
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Figure 23. 


Figure 24. 


Figure 25. 


Figure 26. 


Figure 27. 


Figure 28. 


WIO, and central-eastern tropical Pacific, and the positive 
(negative) OLRAs in the tropical eastern 10 and western tropical 
Pacific, consistent with wet (dry) conditions in the HOA, and 
positive (negative) lOD and EN (LN) conditions. Images created 

using ESRL (2010).57 

Correlations of HOA PR in Oct-Nov with global SST for 1970-2009 
with SST leading by: (a) one month (Sep-Oct), (b) two months 
(Aug-Sep), (c) three months (Jul-Aug), and (d) four months (Jun- 
Jul). The blue box represents the Oct-Nov HOA PR predictand 
region. The red box indicates a region north of New Zealand of 
persistently strong and significant negative correlation between 
SST and the predictand. Note the similarities between the 
correlation patterns in this figure and the SSTA patterns in Figure 

18. Images created using ESRL (2010).59 

Magnitude of correlations between Oct-Nov HOA PR predictand 
and four potential predictors of this predictand, with the predictors 
leading by zero to six months. The four predictors are: NZ SST, 
MEI, Nino3.4, and DMI. The colored text boxes on the right of the 
figure describe the climate variations associated with the last three 
of these potential predictors. Correlations ^ 0.304 indicate 

significance at the 95% level (dashed line).61 

Time series of HOA Oct-Nov PR (mm/day) predictand in Oct-Nov 
(black line) and NZ SST (°C) predictor in (a) Oct-Nov (red line), (b) 
Jul-Aug (purple line), and (c) Apr-May (blue line) during 1970- 
2009. Note the persistence in the SSTs indicated by the similarities 
in the SST time series. Note also the indications of negative 

correlations between the predictand and NZ SSTs.62 

Correlations of HOA PR in Oct-Nov with global SST in Sep-Oct 
based on data from 1970-2009.. These and related correlations 
(Figures 23-24) led us to select SST north of New Zealand (red 
box) as the initial predictor for use in forecasting HOA PR in Oct- 
Nov (our predictand; blue box) at leads of one to six months. 

Image created using ESRL (2010).63 

Average verification metrics by forecast lead time for hindcasts of 
HOA PR in Oct-Nov 1970-2009 using Method 1 (tercile matching) 
and the NZ SST predictor. Averaging is done over all forecast 
categories (AN, NN, and BN PR). Verification metrics consist of 
percent correct (% Corr), false alarm rate (FA Rate), probability of 
detection (POD), and Heidke Skill Score (HSS). Refer to Chapter 
II, Section B.6 for more information on the metrics listed in this 
table. Note the general decline in skill as lead time increases, but 

that HSS remains positive at all lead times.66 

Results from composite analysis forecasts of Oct-Nov HOA PR 
using the NZ SST predictor at a lead time of six months (i.e., using 
NZ SST in Apr-May). (a) Historical frequency distribution of HOA 









Figure 29. 


Figure 30. 


Figure 31. 


Figure 32. 


Figure 33. 


PR probabilities given any SST condition at six months lead. 
Statistically significant results are outlined in black. (b) 
Corresponding probabilistic long-range hindcast of HOA PR in Oct- 
Nov based on AN NZ SST in Apr-May. (c) Corresponding 
probabilistic long-range hindcast of FlOA PR in Oct-Nov based on 

BN NZ SST in Apr-May.68 

Time series by forecast lead time of: (1) correlations between HOA 
PR predictand and NZ SST predictor; and (2) HSS values for tercile 
matching method hindcasts at lead times of zero to six months, and 
using tercile categories based on the most recent 10, 15, and 40 
years, (a) 1995-2009 correlations (green line); (b) 2000-2009 
correlations (red line); (c) 1995-2009 hindcast HSS values (orange 
and light blue lines); and (d) 2000-2009 hindcast HSS values 
(purple and turquoise lines). The bin periods represent the period 
used to form the predictor and predictand AN, NN, and BN 
categories. Note that the correlations shift from negative to 

positive, and the HSS values decrease, as lead time increases.73 

Correlations of HOA PR in Oct-Nov with global SST for 1995-2009 
with SST leading by: (left panel) zero months (Oct-Nov) and (right 
panel) four months (Jun-Jul). The yellow outlined boxes represent 
the multi-region SST predictor regions identified as having the 
highest correlations to the HOA PR predictand at specific lead 
times when using the advanced OCN approach. We used the SST 
predictor east of the Philippines (left panel) was used for 
forecasting at leads of zero to three months. We used the SST 
predictor west of Namibia (right panel) for forecasting at leads of 

four to six months. Images created using ESRL (2010).74 

Composite mean SST (°C) in Apr-May 2009 averaged over: (top 
right) NZ predictor region and (bottom left) MR predictor region. 
The black (red) box represents the NZ (MR) SST predictor regions 
used for analyzing the SSTs at lead times of six months. The area 
averaged SST for Apr-May 2009 was 23.875 (22.958) °C in the NZ 

(MR) predictor regions. Images created using ESRL (2010).78 

Composite mean (left panel) and anomaly (right panel) PR 
(mm/day) in Oct-Nov 2009. The black outlined box represents the 
HOA PR predictand region. Note that AN, NN, and BN PR 
anomalies occurred within the predictand region (right panel); 
however the actual area-averaged PR was 2.836 mm/day, or NN. 

Images created using ESRL (2010).80 

Composite anomalies of PR (mm/day) during Oct-Nov for three 
different years: (left panel) 1997, (middle panel) 1998, and (right 
panel) 2006. The red outlined boxes represent the HOA PR 
predictand region. The HOA experienced extreme AN or BN PR in 
each of these years: AN in 1997 and 2006, and BN in 1998). 
These years were chosen to test the four forecast methods using 
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the non-OCN and OCN-based predictor sets. Images created 

using ESRL (2010).83 

Figure 34. Total verification metrics for hindcasts of HOA PR in Oct-Nov 
averaged over all forecast categories (AN, BN, NN) and all lead 
times of zero to six months. Columns 1-4 (5-8) represent 
verification metrics using the NZ (Multi-Region) SST predictor for 
1970-2009 (1995-2009) with each of the forecast methods listed at 
the bottom right of the figure (e.g., tercile matching, linear 
regression and CAF). For example, column 5 represents the 
verification metrics for Method 1 (tercile matching) using the Multi- 
Region SST predictor. Verification metrics consist of percent 
correct (% Corr), false alarm rate (FA Rate), probability of detection 
(POD), and Fleidke Skill Score (FISS). Refer to Chapter II, Section 
B.6 for more information on these metrics. The bold black line 

highlights the HSSs.88 

Figure 35. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 2 (averaged linear regression) and the NZ 
SST predictor for Oct-Nov 1970-2009 at lead times of zero to six 
months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the periods of relatively 
good model performance before 1984 and after 1993 (see Chapter 
III, Section A.3; Figure 24); and (b) the lower variability in the 
hindcasted PR than the observed PR, especially at longer leads. 

All PR values in mm/day.89 

Figure 36. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 3 (non-averaged linear regression) and the 
NZ SST predictor for Oct-Nov 1970-2009 at lead times of zero to 
six months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the periods of relatively 
good model performance before 1984 and after 1993 (see Chapter 
III, Section A.3; Figure 24); and (b) the lower variability in the 
hindcasted PR than the observed PR, especially at longer leads. 

All PR values in mm/day.90 

Figure 37. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 2 (averaged linear regression) and the MR 
SST predictors for Oct-Nov 1995-2009 at lead times of zero to six 
months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the overall very good model 
performance, especially at shorter leads; and (b) the relatively 
realistic variability, especially at shorter leads. All PR values in 

mm/day.91 

Figure 38. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 3 (non-averaged linear regression) and the 
MR SST predictors for Oct-Nov 1995-2009 at lead times of zero to 
six months. The hindcast PR values are dark red (green) at the 
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Figure 39. 


Figure 40. 


Figure 41. 


Figure 42. 


Figure 43. 


longest (shortest) lead times. Note: (a) the overall very good model 
performance, especially at shorter leads; and (b) the relatively 
realistic variability, especially at shorter leads variability in the 
hindcasted PR than the observed PR, especially at longer leads. 

All PR values in mm/day.92 

Correlations between the actual Oct-Nov HOA PR and the 
hindcasted HOA PR values shown in Figures 35-38 by lead time. 
Forecast method and predictor indicated by colored lines and 

symbols (see key).93 

Conditional composite anomalies of sea surface temperature (°C) 
for the five most extreme BN NZ SST events during Oct-Nov 
1970-2009. Recall that BN NZ SST is associated with EN and 
positive lOD events (Figure 17), and with AN HOA PR during Oct- 
Nov. Red box indicates NZ SST region. Red oval indicates area of 


pronounced negative SSTAs in the tropical Atlantic. Image created 

using ESRL (2010).97 

Conditional composite anomalies of sea surface temperature (°C) 
for the five most extreme EN events (left panel) and LN events 
(right panel) during Oct-Nov 1970-2009. EN and LN events 
identified by use of MEI. Image created using ESRL (2010).98 


(Left panel) 3-D schematic of the anomalous Hadley-Walker 
Circulation (HWC) over the Pacific Ocean corresponding to LN 
events and enhanced convection over the MC. (Right panel) 
Conditional composite anomalies of 850 hPa vector winds (m/s) for 
the five most extreme AN NZ SST events during Oct-Nov 1970- 
2009. Recall that AN NZ SST is associated with LN and negative 
lOD events (Figure 17), and with BN HOA PR during Oct-Nov. 
Bold black line indicates SPCZ. H (L) indicates location of 
anomalous high (low) low-level pressure system. Red arrows 
schematically show the anomalous low-level flow. Image on left 
adapted from Murphree (2009a). Image on right created using 

ESRL (2010).99 

Conditional composite anomalies of 200 hPa GPH (gpm) for the 
five most extreme BN NZ SST events during Oct-Nov 1970-2009. 
Recall that BN NZ SST is associated with EN and positive lOD 
events (Figure 17), and with AN HOA PR during Oct-Nov. The red 
arrows schematically show anomalous wave trains that emanate 
from the maritime continent sector, extend into the extratropcial 
northern and southern hemispheres, and then into the tropics near 
Africa. The arrow directions indicate schematically the direction of 
energy propagation in the anomalous wave trains. The l/l/ef label in 
the upper left indicates that these anomalies are associated with 
above normal precipitation in the HOA. Image created using ESRL 
( 2010 ). 100 
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Figure 44. 


Figure 45. 


Figure 46. 


Figure 47. 


Figure 48. 


Conditional composite anomalies of 200 hPa GPH (gpm) for the 
five most extreme AN NZ SST events during Oct-Nov 1970-2009. 
Recall that AN NZ SST is associated with LN and negative lOD 
events (Figure 17), and with BN FlOA PR during Oct-Nov. The red 
arrows schematically show anomalous wave trains that emanate 
from the maritime continent sector, extend into the extratropical 
northern and southern hemispheres, and then into the tropics near 
Africa. The arrow directions indicate schematically the direction of 
energy propagation in the anomalous wave trains. The Dry label in 
the upper left indicates that these anomalies are associated with 
below normal precipitation in the HOA. Image created using ESRL 

( 2010 ). 101 

Conditional composite anomalies of 850 hPa GPH (gpm) for the 
five most extreme BN NZ SST events during Oct-Nov 1970-2009. 
Recall that BN NZ SST is associated with EN and positive lOD 
events (Figure 17), and with AN HOA PR during Oct-Nov. H (L) 
indicates location of anomalous high (low) low level geopotential 
heights. The red arrows show schematically the anomalous wind 
directions. Note the anomalous convergence into the HOA region. 

Image created using ESRL (2010).103 

Conditional composite anomalies of precipitation rate (PR; mm/day) 
for the five most extreme BN NZ SST events during Oct-Nov 
1970-2009. Recall that BN NZ SST is associated with EN and 
positive lOD events (Figure 17), and with AN HOA PR during Oct- 
Nov. The red box indicates our HOA PR predictand region. Image 

created using ESRL (2010).104 

Conditional composite anomalies of 850 hPa GPH (gpm) for the 
five most extreme BN Namibia SST events during Apr-May 1970- 
2009. Recall that BN Namibia SST is associated with AN HOA PR 
during Oct-Nov. The H indicates an anomalous high-pressure 
region. The direction of the corresponding wind anomalies is 
shown schematically by the red arrow. The cooler SSTs region is 
explained in the main text discussion of Figure 40. Image created 

using ESRL (2010).106 

Conditional composite anomalies of SST (°C) for the five most 
extreme BN Namibia SST events during Apr-May 1970-2009. 
Recall that BN Namibia SST is associated with AN HOA PR during 
Oct-Nov. The negative SSTA region west of Namibia corresponds 
to the Namibia SST predictor region identified as having the highest 
correlations, at leads of four to six months, to the HOA PR 
predictand for 1995-2009 (see Chapter III, Section C.2). The H 
indicates an anomalous lower tropospheric high-pressure region 
and the black arrow indicates the corresponding wind anomaly (see 
Figure 47). Image created using ESRL (2010).107 
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Figure 49. 


Figure 50. 


Figure 51. 


Figure 52. 


Figure 53. 


Correlations of HOA PR in Oct-Nov with global SST for 1970-2009 
(left panel) and 1995-2009 (right panel), with SSTs leading by zero 
months (Oct-Nov). The yellow outlined boxes represent NZ 
predictor region (left panel) and the Philippine SST predictor region 
(right panel). Note the similar patterns for the two periods, but with 
stronger correlations for 1995-2009. Images created using ESRL 

(2010).108 

Composite anomalies of SST (°C) for Oct-Nov of 1995-2009. Note 
the positive SSTAs, especially in the northern hemisphere and in 
the western North and South Pacific. Image created using ESRL 

(2010).109 

Schematic diagram of lower tropospheric and sea surface 

conditions that occur in October-November when: (a) precipitation 
is greater than normal in the FlOA (i.e., during AN FlOA PR events); 
and (b) when SSTs north of New Zealand are cooler than normal 
(i.e., during BN NZ SST events) and. Map background from 

Google Maps (www.maps.google.com) .110 

Schematic diagram of lower tropospheric and sea surface 

conditions that occur in October-November when: (a) precipitation 
is less than normal in the FlOA (i.e., during BN FlOA PR events); 
and (b) when SSTs north of New Zealand are warmer than normal 
(i.e., during AN NZ SST events) and. Map background from 

Google Maps (www.maps.google.com). 111 

Schematic diagram of lower tropospheric and sea surface 

conditions that occur in October-November when precipitation is 
greater than normal in the FlOA (i.e., during AN FlOA PR events). 
Map background from Google Maps (www.maps.google.com).112 
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I. INTRODUCTION 


A. BACKGROUND 

Much of Africa has experienced extreme climate variations in the last 
several decades that have had devastating impacts on the environment and 
societies (Time 2009). The effects of climate change associated with global 
warming have been projected to be especially severe in Africa due in part to on¬ 
going poverty and the lack of strong governance to prepare for and mitigate 
climate change effects (IPCC 2004). Skillful long-range forecasts (LRFs; 
forecasts at lead times of two weeks or longer) of climate variations and climate 
change, and their impacts on specific regions, are critical in managing climate 
related risks, such as the risks of famine, and economic and political instability. 
These LRFs can contribute to the development of early warning systems (EWS), 
such as the Famine Early Warning System Network (FEWS NET; Funk 2009), 
that help military and non-military planners manage the impacts of climate 
variations and climate change. 

In this study, we have investigated the potential for improving LRFs of 
precipitation in the Horn of Africa (HOA). Our over-arching goal is to develop 
long-range forecasting methods that contribute to the development of a more 
holistic approach to dealing with military and non-military challenges in the HOA, 
and to contribute to the welfare of the people of this region—for example, to help 
support early warnings and decisions that save lives and minimize conflicts. 

B. IMPACTS OF CLIMATE ON THE SOCIAL, ECONOMIC, AND 

POLITICAL CONDITIONS IN THE HOA 

The HOA, as defined in this study, comprises six countries (Djibouti, 
Eritrea, Somalia, Ethiopia, Kenya, and Sudan) and covers approximately six 
million square kilometers (ICRC 2004). There are an estimated 260 million 
pastoralists and agro-pastoralists in sub-Saharan Africa, which means 60 to 100 
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percent of the livelihood of most of the regions’ population, and most of the 
region’s gross domestic product comes from agriculture (ICRC 2004). Thus, 
arable land and water (e.g., from rainfall) are extremely important. The 
population is estimated to grow on the order of 100% from 2000 to 2050, leading 
to increasing stress on natural resources and social management structures, and 
increased economic, social, and political tension (ICRC 2004). Competition over 
natural resources is one of the most common causes of conflict in the HOA 
(ICRC 2004). 

Extreme climate events, such as droughts and floods, initiate many 
societal emergencies in the HOA; for example, drought triggers many of these 
emergency events in Ethiopia (ICRC 2004). The effects of these climate related 
emergencies are wide ranging and have multiple direct and indirect cascading 
effects that span the social, economic, political and cultural sectors. Pre-existing 
conditions tend to make the HOA especially vulnerable to the effects of droughts 
and floods. These pre-existing conditions include weak governance, lack of early 
warning and basic social services, poverty, and ongoing recovery from previous 
extreme climate events (UNEP 2009). Extreme climate events often reinforce or 
exacerbate the pre-existing conditions (UNEP 2009). Some of the harsh 
negative effects that are the direct or indirect results of droughts and floods are 
(e.g., ICRC 2004; UNEP 2009; Time 2009; Ogallo 2010): 

(1) Shortages of food and water 

(2) Environmental degradation 

(3) Large losses to hydroelectric power generation 

(4) Agriculture, wildlife and livestock devastation 

(5) Famine and starvation 

(6) Increased health risks and incidences of diseases 

(7) Internally displaced persons (IDPs), refugees, and forced migration 

(8) Large economic losses, trade bans, and border restrictions 
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(9) Economic and political mismanagement; strongman governments 

(10) Border conflicts and extremism 

This list of effects is not all-inclusive and does not capture the interplay 
between effects that can greatly amplify the environmental, social, economic, and 
political impacts of climate variations. The hazards of extreme climate events are 
highlighted by several recent cases. The prolonged droughts in the Darfur region 
of Sudan during recent decades that have greatly contributed to conflicts over 
control of scarce resources, such as fertile land and water, that have led in turn 
to over 300,000 deaths and two million DPs since 2003 (UNEP 2009). During 
1992-1993, famine related to drought helped create political instability that let the 
U.S. to send troops into Somalia; the resulting conflicts led to the deaths of 1,000 
Somalis and 18 Americans (Baker 2010). In the past decade, approximately 23 
million people in seven East African countries have been fed by aid agencies 
after a decade of poor rains have led to food and water shortages (CNN 2009; 
Time 2009). In 2006, the UN said it expected climate change related to global 
warming to have the greatest societal impacts in Africa, where nearly 90 million 
Africans are especially vulnerable to the impacts of extreme weather and climate 
events, such as droughts and floods (Time 2009). 

C. U.S. AND INTERNATIONAL POLICY TOWARD HOA 

The U.S., UN, aid agencies, and other nations and organizations have 
begun to account for climate factors in their attempts to manage environmental, 
societal, economic, and political problems in the HOA and other regions (e.g., 
UNEP 2009). The long-standing band-aid and reactive responses climate 
related problems are being replaced by early warning systems (EWSs) designed 
to integrate environmental analyses and forecasts into societal, economic, and 
political planning, and conflict prevention and peace building efforts. The UN has 
strongly recommended increased development of EWSs and early action in 
countries that are vulnerable to conflicts over natural resources and 
environmental issues (UNEP 2009). U.S. policy toward such countries is 
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evolving from primarily a crisis intervention approach to an approach based on 
identifying and responding to the factors that cause the crises (e.g., Baker 2010; 
U.S Africa Command 2010). 

The U.S. recognizes the importance of early action in the HOA because 
the region is prone to extremist groups and predatory elites that can thrive largely 
with impunity, and that often exacerbate the suffering of the populations under 
their control (Baker 2010). The Army Field Manual 3-07 (USAFM 2008), states 
the greatest threat to U.S. national security will not come from emerging 
ambitious states but from nations unable or unwilling to meet the basic needs 
and aspirations of their people (Baker 2010). The Failed States Index of 2009 
listed five out of the six HOA countries in alert status; Somalia and Sudan being 
number one and three, respectively (FSI 2009). Extreme climate events can 
escalate conflicts and generate others, especially in failed or failing states (UNEP 
2009). Thus, investments in EWSs for conflict prevention and peacebuilding, in 
addition to humanitarian assistance, are especially important in nations that are 
in, or vulnerable to, failed state status. Better long-range forecasts of 
environmental conditions (e.g., of precipitation in the HOA) are vital to the 
operation of EWSs and other aspects of a holistic approach to dealing with 
climate related societal, economic, and political problems. 

D. CLIMATE VARIATIONS AND TELECONNECTIONS IN THE HOA 

Before attempting to develop methods for long-range forecasting of 
climate variations in the HOA, it is important to understand what the major 
variations are and what regional and global scale dynamical processes help 
create these variations. Numerous prior studies have identified intraseasonal to 
interdecadal variations that tend to occur in the HOA and nearby regions (e.g., 
Nicholson and Entekhabi 1987; Farmer 1988; Ogallo 1988; Hutchinson 1992; 
Hastenrath et al. 1993, 2004; Mutai et al. 1998; Kabanda and Jury 1999; Mutai 
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and Ward 2000; Poccard et al. 2000; Camberlin et al. 2001; Black et al. 2002; 
Clark et al. 2002; Behera et al. 2005; LaJoie 2006; Vorhees 2006; Montgomery 
2007). 

The focus of many of these studies has been on interannual climate 
variations in precipitation, in large part because these variations can be very 
pronounced and have serious impacts on the environment and the human 
population. These precipitation variations are generally expressed as extremes 
in the precipitation that occurs during the two main rainy seasons in east African 
and the HOA region—the long rains during Mar-May and the short rains during 
Oct-Dec (e.g., Mutai and Ward 2000). The precipitation variations tend to be 
associated with teleconnections between HOA precipitation and sea surface 
temperature (SST) in other regions of the tropics and subtropics (e.g., Ogallo 
1988; Hutchinson 1992). Some of the most notable teleconnections are those 
associated with the SST changes that occur in the tropical Pacific and Indian 
basins during El Nino—La Nina (ENLN) and Indian Ocean Dipole (lOD; also 
known as Indian Ocean Zonal Mode (lOZM); Saji et al. 1999) events. Numerous 
studies have linked variations in Oct-Nov precipitation in the HOA with ENLN, 
lOD, and, to a lesser extent, with SST variations in the Atlantic Ocean (e.g., 
Nicholson and Entekhabi 1987; Farmer 1988; Poccard et al. 2000; Camberlin et 
al. 2001; Black et al. 2002; Clark et al. 2002; Behera et al. 2005; LaJoie 2006; 
Vorhees 2006). 

Much of the impact of remote SSTs on HOA precipitation occurs via 
fluctuations of the circulation over the tropical and subtropical Indian Ocean (10). 
As SSTs decrease (increase) near the maritime continent during EN (LN) events, 
the zonally oriented Walker Circulation over the 10 tends to weaken (strengthen), 
which results in: (a) less (more) subsidence over east Africa and the HOA; (b) 
more (less) subsidence over the maritime continent; (c) weaker (stronger) low 
level westerlies over the tropical 10; (d) warmer (cooler) SSTs off east Africa and 
the HOA; (e) cooler (warmer) SSTs in the eastern tropical 10; (f) more (less) 
favorable conditions for precipitation over east Africa, the HOA, and the western 
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topical 10; and (g) less (more) favorable conditions for precipitation over the 
maritime continent and the eastern tropical 10 (e.g., Ogallo 1988; Mutai and 
Ward 2000; Black et al. 2002; Clark and Cole 2002; LaJoie 2006; Vorhees 2006). 
A similar set of variations tends to occur during the positive (negative) phases of 
the lOD (Saji et al. 1999; Behera et al. 2005; Vorhees 2006; Twigg 2008). The 
lOD and ENLN are well correlated, with the positive (negative) lOD phase 
tending to occur during EN (LN) events (Twigg 2008) ENLN and lOD events 
tend to reach their maximum amplitudes during September-January. Thus, they 
tend to have relatively large impacts on the short rains during Oct-Dec ((cf. 
Hastenrath et al. 1993; LaJoie 2006; Vorhees 2006). Figure 1 provides a 
schematic illustration of the typical relationships between HOA precipitation and 
EN and positive lOD events. 

The typical relationships between HOA precipitation variations and ENLN 
and lOD events can break down if ENLN and lOD do not follow their usual 
relationships (e.g., EN with positive lOD; LN with negative lOD), or if they peak at 
different times than is typical, or if they occur simultaneously with other climate 
variations (e.g., Madden-Julian Oscillation; MJO) (cf. LaJoie 2006; Vorhees 
2006). 

Public awareness to these climate variations is increasing as mainstream 
media highlight the associated environmental, societal, and economic impacts of 
the variations. Figure 1 is an example of how the news media is informing the 
general public about ENLN, lOD, and their global scale impacts. 


6 




Figure 1. Schematic illustration of anomalies in the Indian and Pacific basin 
Walker circulations, and in low-level atmospheric pressures and 
sea-surface temperatures, during EN and positive lOD phase 
events. These anomalous conditions tend to contribute to above 
normal precipitation in the HOA. These anomalies tend to be most 
pronounced during August-February, and thus tend to have their 
greatest impact on the HOA short rains period in October- 
December. Figure from Hunt (2009). 


Teleconnections between the east Africa and HOA short rains and 
conditions in the Atlantic Ocean are also important, but have not been as 
extensively studied as teleconnections between the Pacific and lO basins. 
Nicholson and Entekhabi (1987) found that below (above) normal SSTs in the 
equatorial Atlantic off the coast of West Africa tended to be associated with 
below (above) normal precipitation in the east Africa short rains period. 
However, Mutai and Ward (2000) found that above (below) normal SSTs in the 
equatorial and tropical Atlantic off the coast of Africa were associated with above 
(below) normal east African short rains. 
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E. NON-DOD OPERATIONS IN THE HOA AND LONG-RANGE 

FORECASTING 

The regional and international community has instituted EWSs and other 
preventative systems for much of Africa and the HOA in particular. EWSs are 
vital to the HOA because the early identification of populations at risk can enable 
the timely focused actions needed to avert widespread hunger, destitution, 
famine, and violent conflict. Several EWSs use the results from climate 
monitoring and long-range forecasting to develop early warnings of food and 
agricultural insecurity, drought, flooding, and famine (Verdin et al. 2005). 

The Famine Early Warning System Network (FEWS NET), an activity of 
the U.S. Agency for International Development (USAID), is a multi-disciplinary 
project that collects, analyzes, and distributes national and regional information 
to decision makers about potential or current climate, famine, or socio-economic 
hazards (FEWS NET 2010). This allows decision makers to authorize timely 
measures to prevent food insecurity and other societal problems related to 
climate factors (e.g., drough, flooding, high temperatures). The U.S. National 
Aeronautics and Space Administration (NASA), the U.S. Geological Survey 
(USGS), U.S. Department of Agriculture (USDA), and NOAA provide satellite 
based analyses of recent HOA temperature, precipitation, and land surface 
vegetation, plus other scientific and technical support, to FEWS NET (FEWS 
NET 2010). These analyses and related long-range forecasts of climate 
conditions (e.g., surface air temperature, precipitation) are used to develop food 
insecurity outlooks at lead times of several months or less. In particular, the 
environmental analyses and forecasts are combined with in-country analyses of 
prices, grain stores, political conditions, and agriculture inputs to provide 
assessments of present and potential food shortages (Funk 2009). Figure 2 
outlines the process used to generate the climate forecasts. Figure 3 outlines 
the process used to develop a food security outlook product based on climate 
analyses and forecasts, and on inputs from experts on climate, agriculture. 
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health, economics, disaster management, and many other fields. Figure 4 
provides an example of a food security outlook product from FEWS NET. This 
type of outlook is intended for use by public and private organizations in planning 
mitigation response activities. 




Precipitation 

Anomaly 

Forecasts 


Figure 2. Outline of process used to generate some of the climate forecasts 
used by FEWS NET. Observations (left box) are used to identify 
statistical relationships between climate system variables (middle 
box; e.g., conditions in the Indian Ocean) and food related variables 
(e.g., precipitation in southern Africa). These relationships are 
used to build predictive statistical models that generate forecasts 
for the food security related variables (right box). Figure from Funk 
(2009). 



Figure 3. Outline of contingency planning and response process used by 
FEWS NET. Climate forecasts are used as inputs to the 
development of potential food insecurity scenarios that are then 
used to develop contingency plans and food security outlooks (e.g.. 
Figure 4). Figure from Funk (2009). 
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Figure 4. Example of a food security outlook product generated by FEWS 
NET. This outlook is for conditions in eastern central Africa and the 
FlOA in April-June 2010. Climate forecasts are an essential 
element in the generation of such outlooks. These outlooks are 
used by agricultural, emergency management, and other planners 
to mitigate the food insecurity events and their effects. Figure from 
U.S. Aid, FEWSNET, May 2010; available online at 
http://www.fews.net/Pages/region.aspx?gb=r2&l=en; accessed May 
2010 . 


Other organizations that incorporate climate forecasts into their EWSs 
include the Intergovernmental Authority on Development (IGAD) Climate 
Prediction and Applications Center (ICPAC; IGAD 2010). The Drought 
Monitoring Center (DMC) in Nairobi, Kenya, a specialized institution of the 
ICPAC, is responsible for climate monitoring, prediction, early warning, and 
applications for the reduction of climate-related risks in the greater HOA (WMO 
2010a, b). ICPAC provides regular climate advisories, including 10-day, monthly, 
and seasonal climate bulletins as well as early warning information on evolving 
climate extremes and associated impacts. Prior to each rainy season a greater 
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HOA climate outlook forum is held that brings together international, regional, 
and national experts in climate, health, agriculture, economics, and disaster 
management. This forum examines prevailing and expected climate variations 
that affect weather and climate in the HOA. Using atmosphere-ocean dynamical 
models, statistical models, and expert interpretation, a consensus regional 
seasonal forecast is produced in the form of probability distributions of above 
normal (AN; upper tercile), near normal (NN; middle tercile), and below normal 
(BN; lower tercile) precipitation for the greater HOA. The contributors to the 
regional consensus climate outlook include the ICPAC regional meteorological 
offices. International Research Institute for Climate and Society (IRI), FEWS 
NET, United Kingdom Meteorological Office (UKMO), World Meteorological 
Organization (WMO), National Centers for Environmental Prediction (NCEP) / 
Climate Prediction Center (CPC), and the European Centre for Medium-range 
Weather Forecasts (ECMWF) (ICPAC 2010). Figure 5 is an example of the 
ICPAC precipitation outlook for the greater HOA for Mar-May 2010. 
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Longitude (*E) 

Figure 5. Example of a seasonal precipitation outlook from ICPAC. This 
outlook shows the predicted probability of above normal (AN), near 
normal (NN), and below normal (BN) precipitation in the greater 
Horn of Africa during March-May 2010. The yellow, green, and 
gray shading represents areas expected to experience similar AN, 
NN, and BN precipitation conditions. Yellow: relatively high 

probability of NN or BN precipitation. Green: relatively high 

probability of NN or AN conditions. Gray: equal chances of AN, 
NN, and BN. The vertical boxes with three numbers show the 
probabilities for each category, with AN at the top of the box, NN in 
the middle, and BN at the bottom. Note that in each colored region, 
the three probabilities sum to 100%. Example: in zone III the 
probabilities are 35% AN; 40% NN; and 25% BN. Figure from 
ICPAC; available online at www.icpac.net/Forecasts; accessed 
March 2010. 
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The Kenya Meteorological Department (KMD) produces seasonal 
precipitation forecasts for Kenya of precipitation probabilities, and of the dates of 
the onset and cessation of precipitation. These forecasts are based mainly on 
empirical statistical models for which the expected evolution of global SST 
anomalies and the Southern Oscillation Index (SOI) are the main predictors 
(KMD 2009). Figure 6 shows examples of KMD seasonal forecasts for Kenyan 
precipitation for March-May 2010. 
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Figure 6. Examples of seasonal precipitation forecasts from the Kenya 
Meteorological Department for March-May 2010. March-May in 
central east Africa is normally a period of relatively high 
precipitation that is often referred to as the long rains season. Top 
panel: forecast of areas of above and below normal precipitation. 
Bottom left panel: forecast of onset dates for long rains season. 
Bottom right panel: forecast of cessation dates for long rains 
season. Figure from Kenya Meteorological Department; available 
online at www.meteo.qo.ke/ ; accessed March 2010. 
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F. NON-DOD OPERATIONS IN THE HOA AND LONG-RANGE 

FORECASTING 

DoD has recently been restructured in order to address African peace and 
stability issues, including the creation of U.S. Africa Command, or AFRICOM. 
AFRICOM is directly responsible to the U.S. Secretary of Defense and is the only 
U.S. combatant commands devoted solely to Africa. Part of the AFRICOM 
mission is to enable DoD and other elements of the U.S. government to work in 
concert and with partners to achieve a more stable environment in which political 
and economic growth can take place (U.S. Africa Command 2010). Unlike 
traditional combatant commands, AFRICOM’s focus is on war prevention rather 
than on war-fighting, including working with over 53 African nations and 
organizations to build regional stability and crises response capacity in support of 
U.S. government efforts in Africa. Among President Obama’s main priorities for 
Africa are: (a) helping to develop strong and sustainable democratic 
governments; (b) strengthening public health; and (c) promoting peace and 
security (U.S. Africa Command 2010). Some of the AFRICOM activities include 
health services, humanitarian aid and support, and training of the African military 
and security forces—all in an effort to support the African people in peacefully 
building a foundation on which to grow economically and politically. Figure 7 
illustrates some of these activities, all of which benefit to some degree from 
accurate long-range forecasts of precipitation and other climate system variables. 
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Figure 7. Examples of U.S. Africa Command (AFRICOM) activities that 
contribute to U.S. government efforts in Africa. Accurate long- 
range climate forecasts have the potential to improve planning by 
AFRICOM and its African government partners and non¬ 
governmental organization (NGO) partners. Figure from U.S. Africa 
Command (2010). 


If AFRICOM is to be successful, it must account for the many climate 
factors that affect Africa, including climate variations and climate change. To do 
so, it will have to collaborate with operational weather and climate organizations 
that have the capacity to provide climate analyses and long-range predictions of 
climate factors that affect operations in Africa (Montgomery 2008). Presently, 
AFRICOM relies mainly and indirectly on civilian climate organizations (e.g., 
NOAA, UKMO). Flowever, the US Air Force (USAF) 14th Weather Squadron (14 
WS) has recently started producing long-range forecasts (LRFs) for the FlOA that 
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support military operations in the region. The 14 WS is also investigating the 
possibility of producing similar forecasts for the pan-sub-Saharan African region 
(e.g., the Sahel). 

The 14th Weather Squadron (14 WS), formerly known as the Air Force 
Combat Climatology Center, located in Asheville, NC is the primary operational 
climate support center for the USAF and much of the rest of DoD (14*^ Weather 
Squadron 2010). Its mission is to receive, store, and quality control worldwide 
weather data to create climatological products to strengthen the combat 
capability of U.S. warfighters. Customers include all the military services, the 
combatant commands (e.g., AFRICOM), the USAF MAJCOMS, U.S. intelligence 
agencies, the White House, and many more. Military planners intending to 
exploit the environment for operations worldwide, including the HOA region, use 
the information and services provided by the 14 WS. The information in this 
section on the 14 WS LRFs is based on personal communication with the 14 WS 
personnel and examination of sample LRFs from the 14 WS. 

In January 2010, the 14 WS started producing narrative-only LRFs (lead 
times of one to six months) of HOA surface air temperature, precipitation, and 
cloud ceilings. The forecasts include discussions of: (a) the state of El Nino - La 
Nina (ENLN); (b) expected departures from normal conditions; (c) forecast 
summary conditions; and (d) annual climatologies for each country within the 
HOA. Figure 8 shows the region for which the 14 WS produces its HOA LRFs. 
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Figure 8. The African region (outlined in yellow) for which the 14th Weather 
Squadron (14 WS) produces its LRFs for the FIOA. The 14 WS 
produces for this region narrative-only long-range forecasts (lead 
times of one to six months) of surface air temperature, precipitation, 
and cloud ceilings. These LRFs are based on composite analysis 
methods employed by the 14 WS, and on LRFs produced by 
civilian operational climate centers (e.g., CPC, IRI, UKMO). The 
locations marked with red dots are ones from which in situ 
observations are collected and for which terminal aerodrome 
forecasts (TAFs) are issued. Figure from personal communication 
with 14 WS. 


To produce these LRFs, the 14 WS generates composite analysis 
forecasts (CAFs; discussed in Chapter II) that are then synthesized with a suite 
of climate analysis and forecast products from CPC, IRI, UKMO, Australian 
Bureau of Meteorology (BOM), Beijing Climate Center, Flydrometeorological 
Center of Russia, University College London (UCL) Dept of Space and Climate 
Physics, USDA, and Air Force Weather Agency (AFWA) dust transport products. 
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The 14 WS uses the composite analysis forecast (CAF) process to develop LRFs 
based on using Nino 3.4 sea surface temperature (SST) anomalies (SSTAs) as a 
predictor. University of Delaware (or NCEP reanalysis) precipitation data is used 
to develop the HOA precipitation predictand data set. The resulting CAF LRFs 
predict the probability of AN, NN, and BN conditions for the FlOA for a given 
month. Monthly and seasonal prediction products from outside sources (e.g., 
CPC, IRI, UKMO, BOM) are evaluated and synthesized with the 14 WS CAF LRF 
to produce the final narrative-only LRFs of surface air temperature, precipitation, 
and cloud ceilings. The narrative also discusses the potential for droughts, 
floods, and dust storms, and expected snow pack conditions. The LRF is 
updated monthly or on a specific needs basis valid for 1-180 days from 
issuance. 

The 14 WS uses NCEP reanalysis composites to verify their forecasts. 
The composites are subjectively analyzed and categorized into AN, NN, or BN 
terciles for several sub-regions of the HOA and compared to the forecasts made 
90, 60, and 30 days prior. For example, a 90-day forecast made in January valid 
for April is compared to the April reanalysis composite. The subjective tercile 
categorizations from both the forecast and observed data are manually fed into a 
program that computes the Heidke skill and Brier skill scores, the primary 
verification metrics for the LRFs. 

Long-range forecasts, narratives and other products and services from the 
14 WS can be found at https://notus2.afccc.af.mil/SCIS/ . The user will need a 
password or U.S. government common access card to view the suite of products. 

To the best of our knowledge, no other military service provides LRFs for 
the HOA at the time of this research. 
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G. 


MOTIVATIONS FOR AND OUTLINE OF THIS STUDY 


1. Scientific Motivations 

We had several scientific and operational motivations for conducting this 
research. Scientifically, very few climate studies focus on the HOA as defined in 
this study; instead, focusing further south in central East Africa (especially 
Kenya, Tanzania, Uganda, Zimbabwe; e.g.. Funk 2009). Vorhees (2006) and 
LaJoie (2006) conducted analyses of the patterns, processes, and 
teleconnections associated with precipitation and temperature anomalies in the 
HOA and nearby regions. They identified physical linkages between precipitation 
variability and global scale climate variations, especially ENLN, the lOZM, and 
the MJO. They both noted the potential for LRFs of HOA surface air temperature 
and precipitation, but did not conduct any forecasting research. Relatively little 
research has been done on long range prediction for the HOA, and most of these 
studies have focused on seasonally averaged forecasts at lead times of three 
months or less. Many of the LRF studies for the HOA and East Africa have used 
SOI, ENLN, lOZM and other forecasted and well-known climate indices as 
predictors. Heidt (2009) developed deterministic and probabilistic LRFs for the 
western North Pacific that were based on SST predictors that were specifically 
selected for the predictand variables, times, and locations. Our study builds on 
the climate analysis results of Vorhees (2006) and LaJoie (2006), plus the 
forecasting methodology used by Heidt (2009), to identify and apply predictors 
that are specific for LRFs of HOA precipitation. We also develop and test 
methods for skillful forecasting at leads of up to six months. 

2. Operational Motivations 

With the recent advent of long-range forecast support for the HOA by Air 
Force Weather (AFW), the 14 WS has requested the assistance of NPS in 
improving their long-range forecasting for the HOA, using quantitative, objective 
methods with demonstrated skill. At present, the 14 WS LRFs for the HOA are 
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narrative-only, lengthy, generalized, have lower skill than desired, and are not 
well suited for use in warfighter planning. Part of our objective in this study was 
to develop and test methods for the 14 WS to produce more reliable and 
actionable LRFs. We also intended to improve the climate forecasts used by 
governmental and non-governmental planners and governmental organizations 
in preparing for future environmental, agricultural, societal, and economic 
impacts, and in mitigating social and political instability in the HOA. The vast 
majority of the people living in the HOA depend heavily on precipitation for their 
livelihoods, so better LRFs of precipitation, especially during the relatively 
unpredictable short rain period, would be very beneficial to them. 

3. Research Questions 

This study explored the viability of using advanced climate datasets and 
methods to skillfully forecast atmospheric conditions at long lead times, in order 
to provide DoD and non-DoD organizations with cost-effective information to 
significantly enhance their planning processes. We focused primarily on 
investigating the following questions: 

(1) What atmospheric variables in the HOA are both operationally 
significant and predictable at long lead times? 

(2) What climate system variables are the most viable predictors of climate 
variations in the HOA, and can these variables be used to skillfully predict 
atmospheric conditions in the HOA at long lead times? 

(3) What are the best LRF methods to use at all lead times in the HOA? 

(4) What are the best formats for presenting LRFs to planners preparing 
for operations in the HOA? 

4. Study Region 

Our general study region was the broad HOA Sudan, Ethiopia, Eritrea, 
Djibouti, Somalia, and Kenya and the adjacent waters (Figure 9). 
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Figure 9. The broad Horn of Africa (HOA) region. The continental scale view 
(left panel) shows the HOA region (red outlined box) in relation to 
all of Africa. The HOA focused view (right panel) shows the 
general HOA region for our study. The darker brown shading 
indicates the countries that are often considered part of the broad 
HOA region. Figure from World Health Organization (WHO); 
available online at: www.who.int/—/hoafrica/en/index.html : 
accessed May 2010. 


We chose to focus on this region for several reasons. This region is an 
important area of operations for DoD and non-DoD organizations, especially due 
to political instability, and terrorist and piracy activities, in the region. There are 
extensive interactions of governmental and non-governmental organizations in 
the region (e.g., between Ethiopia, Somalia, and Kenya). There are a number of 
broad socio-economic similarities for the countries in this region (e.g., similar 
agricultural practices, economies, religions). There are also many environmental 
similarities (e.g., similar geography, weather, climate, land surface vegetation). 
In particular, the region tends to experience similar seasonal cycles and 
variations in climate (e.g., monsoonal cycles, rainy seasons, intraseasonal to 
interannual variations in temperature and precipitation; Figure 11). 
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For most of the HOA, precipitation tends to be focused in two high 
precipitation periods: the March-May long rains period and the October- 
December short rains period (Figures 10, 11). These rainy periods occur as the 
intertropical convergence zone (ITCZ) passes overhead as it shifts northward 
(southward) over eastern Africa during the northern spring and summer (fall and 
winter). Figure 10 shows the precipitation patterns associated with the 
southward shift of the ITCZ during the northern fall. 



1 i J * 3 O / tS 


Figure 10. Long-term mean (LTM) composite of surface precipitation rate (PR; 

mm/day) during October-November (Oct-Nov) for 1970-2009. 
Oct-Nov tends to be a period of relatively high precipitation for the 
general FlOA region (see Figures 9 and 11). HOA precipitation in 
this period is: (1) lower than in the western and central tropical 
African and tropical Indian Ocean regions; and (2) higher in these 
months that in the tropical and subtropical regions to the north and 
south of the HOA. Image created using ESRL (2010). 
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5. 


Period 


The focus period for this study was October-November of the 40-year 
period 1970-2009. This period was chosen based on the availability and quality 
of NCEP reanalysis and other datasets. These datasets are available for earlier 
years, but we chose to maximize the positive impacts of satellite data on the 
reanalysis data by excluding years prior to 1970. The 1970-2009 period is long 
enough to provide relatively good depictions of intraseasonal to interannual 
climate variations, such as the Madden-Julian Oscillation (MJO), ENLN, and the 
Indian Ocean Zonal Mode (lOZM) (e.g., Vorhees 2006), as well as some 
representation of decadal variations. We focused on October-November 
because this short rains period provides much of the annual precipitation for the 
HOA and is the most rainy period for the HOA (e.g., Hastenrath et al. 1993; 
Ogallo 1988; Behera et al. 2005), as shown in Figure 11. 



Figure 11. Monthly precipitation rate (PR; mm/day; black bars) and standard 
deviation of the precipitation rate (SD; mm/day; white bars) in the 
east African region (5°S-5°N, 35°-46°E) from in situ observations. 
Two rainy seasons occur in east Africa, including the HOA as 
defined in this study: (1) the March to May long rains period; and 
(2) the October to December short rains period. The PR during 
Oct-Nov tends to be relatively high and has a relatively high SD. 
Thus, the Oct-Nov PR is an especially important, but difficult, 
variable to be able to predict at long lead times. The Oct-Nov PR 
in the HOA was the primary target, or predictand, for the LRFs we 
developed and tested in this study. Figure from Behera et al. 2005. 
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6. Thesis Organization 


Chapter II provides an overview of the datasets, and analysis and LRF 
methods, used in this study. Chapter III outlines the results of each LRF method 
when applied to various periods and predictors, and concludes with a LRF of 
precipitation rate (PR) in the HOA during October-November 2010. Chapter IV 
provides a summary of our results and conclusions, and offers suggestions for 
future research. 
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II. DATA AND METHODS 


A. DATASETS AND SOURCES 

1. NCEP/NCAR Atmospheric Reanalysis Data 

The main climate dataset we used in this study is the NCEP/NCAR 
reanalysis dataset (R1: Kalnay et al. 1996; Kistler et al. 2001). This dataset 
provides a global retrospective analysis (i.e., a reanalysis) of atmospheric and 
surface climate system variables from 1948 to present. Observational data from 
in situ, airborne, and satellite sensors are quality controlled and assimilated with 
a data assimilation system that is kept unchanged over the reanalysis period. 
We used reanalysis data at the standard temporal resolution of six hours and 
horizontal resolution of 2.5°, and at standard tropospheric and stratospheric 
levels. Our primary variables of interest were air temperature, sea surface 
temperatures (SST), geopotential heights (GPH), winds, and precipitation rate 
(PR). 

The strengths of the R1 dataset are its uniform global coverage; 
accessibility; and ability to capture low frequency climate variations. However, 
the 2.5°spatial resolution limits our ability to represent mesoscale features, such 
as those associated with the HOA’s complex topography. 

Our study period of 1970 to 2009 was chosen for the reasons given in 
Chapter I, Section D.2, and because prior studies have shown the R1 dataset is 
problematic and may be unreliable for tropical Africa prior to 1970 (Pocard et al. 
2000; Camberlin et al. 2001). 

NCEP released a major new coupled atmosphere-ocean-land-ice 
reanalysis dataset, the Climate Forecast System Reanalysis (CFSR), in January 
2010 that includes major improvements to the data assimilation and underlying 
dynamical model. CFSR spans the period January 1979-present. Some of the 
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improvements over prior reanalyses are a) coupling to the ocean during the 
generation of the six-hour guess field, b) an interactive sea-ice model, c) 
assimilation of satellite radiances for the entire period and d) higher spatial (0.5° 
X 0.5° in the atmosphere) and temporal (hourly) resolution (Saha et al. 2010). 
CFSR became available too late in our study to allow us to use it. However, 
CFSR should be considered for use in similar future studies. 

2. Multivariate ENSO Index (MEI) 

We analyzed the MEI as a possible predictor given the significant links 
between the HOA and ENLN (e.g.. Black et al. 2002; Behera et al. 2005; LaJoie 
2006; Vorhees 2006). The MEI combines six atmospheric and oceanic variables 
into one index to characterize and monitor current and past states of ENLN in the 
tropical Pacific (Wolter and Timlin 1993; Wolter and Timlin 1998). Because the 
MEI integrates multiple variables into a single bi-monthly index, the MEI is a 
more integrated and stable index of ENLN and a potentially more useful predictor 
than other ENLN indices that only monitor one variable (e.g., SST from the Nino 
3.4 index or sea level pressure (SLP from the Southern Oscillation Index (SOI)). 
A positive (negative) value of MEI represents an EN (LN). 

3. Nino 3.4 Index (Nino3.4) 

We also investigated as a potential predictor of HOA PR another widely 
used index of ENLN (Trenberth 1996), the Nino 3.4 index (Nino3.4) based on 
area-averaged SST anomalies (SSTAs) in the region 5°N - 5°S, 120° - 170°W. 
This index is defined as the three-month average of the SST anomalies, with 
index values greater than or equal to +0.5 °C (less than or equal to -0.5 °C) 
indicating EN (LN) conditions.. We investigated the use of the Nino3.4 as a 
predictor of HOA PR because 1) the 14 WS uses this index as its primary 
predictor for the HOA, and 2) long-range probabilistic forecasts of this index are 
available from CPC and the International Research Institute for Climate and 
Society (IRI). 
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4. 


Indian Ocean Dipole Mode Index (DMI) 


We also investigated as a potential predictor of HOA PR an index of the 
Indian Ocean Dipole Mode, or Indian Ocean Zonal Mode (Saji et al. 1999). For 
this, we used the Dipole Mode index (DMI) based on SSTAs in the Indian Ocean. 
The DMI is based on the SSTA difference between the western equatorial Indian 
Ocean (50°E-70°E and 10°S-10°N) and the southeastern equatorial Indian 
Ocean (90°E-110°E and 10°S-0°N). A positive (negative) DMI represents a 
positive (negative) phase lOD (Saji et al. 1999) and is associated with above 
(below) normal PR in the HOA during October-November (e.g., Behera et al. 
2005; Vorhees 2006). The DMI is also forecasted, with long-range forecasts 
available from the Bureau of Meteorology (BOM) in Australia and the Japanese 
Agency for Marine-Earth Science and Technology (JAMSTEC), from which all 
DMI data for this study were obtained. 

B. ANALYSIS AND FORECASTING METHODS 

Much of the analysis and forecast methodology for this study is similar to 
that used by Heidt (2009) in her study of LRFs of ocean conditions in the western 
North Pacific. 

1. Predictand Selection 

The predictand, or forecast target, for our study was area-averaged 
precipitation rate (PR) in October-November for a selected region within the 
HOA at lead times of zero to six months. The major reasons for selecting this 
predictand are given in Chapter I, Section D. A primary motivation for using 
area-averaged predictands is that they are simpler and larger forecast targets 
than point forecast targets, which tends to: (a) make the development of the 
forecast method simpler; (b) increase predictability at long lead times (e.g., leads 
of several months or seasons); and (c) simplify forecast verification (cf. van den 
Dool 2007). The disadvantage of using a spatially large predictand is that the 
forecast applies uniformly to the whole predictand region and does not identify 
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variability within the region. However, this disadvantage can be mitigated by 
selecting a predictand region that tends to be spatially coherent (e.g., a region for 
which the temporal variations of the predictand variable tend to be spatially 
uniform). Also, this disadvantage may not be as important to end users of the 
forecasts as they are to the forecasters. This tends to be the case for operational 
planners who, at long lead times, only need LRFs for a general region. The need 
for such threshold forecasts is common for DoD planners, who often only need to 
develop spatially broad plans at long leads, with the spatial details being planned 
at shorter leads. Finally, if LRFs of more specific forecast regions do not have 
skill, then spatially broad LRFs may be the only viable option. 

The major factors involved in our selection of the specific predictand 
region within the HOA were: 

(1) The spatial patterns of environmental variability in and near the region 

(2) Significant correlations between PR in the region and potential 
predictors at lags of 0-6 months 

(3) The operational importance of the region 

(4) The need for long-range forecast support in the region 

(5) The overlap with the HOA region defined by the 14 WS 

To select the specific region within the HOA, we began by comparing the 
LTM October-November (Oct-Nov) PR composite (Figure 10) with the PR 
anomaly for October-November of each year during 1970-2009 (not shown). 
These comparisons allowed us to identify a region within which: (a) the LTM PR 
and the PR anomalies were relatively homogeneous, especially for the strongest 
positive and negative PR anomalies. We first conducted these comparisons for 
the broad HOA region bounded by 5S-20N, 30-55E, which encompasses all of 
Ethiopia, Eritrea, Djibouti, Yemen, Somalia, and Kenya. This region also 
includes most of the HOA region defined by the 14 WS (Figure 8). Based on 
these initial comparisons, we subdivided this broad HOA region into 20 sub- 
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regions for which we constructed PR time series that we then correlated with 
each other and with potential predictors (e.g., MEI, Nino3.4, DMI, global SST). 
From this process, we identified the sub-regions for which: (a) the PR and PR 
anomaly patterns were large, uniform, and spatially coherent; and (b) the PR 
correlations with potential predictors were large and statistically significant. 

Based on this process, we selected as our predictand region, the 
rectangular area bounded by 5S-10N, 30-55E. This predictand includes almost 
all of what is generally considered to be the HOA, but also includes some 
adjacent land and ocean areas. We did not include the northern portion of the 
broad HOA region (e.g., northern Ethiopia, Eritrea, Djibouti, and most of Sudan) 
in our predictand region because it was too dissimilar from the HOA region 
further to the south, primarily in terms of the signs, magnitudes, and spatial and 
temporal patterns of the PR anomalies. 

2. Composite Anomalies, Correlations, and Teleconnections 

We used the predictand region PR times series for Oct-Nov of 1970-2009 
to identify: (a) the above normal (AN), near normal (NN), and below normal (BN) 
terciles (i.e., the upper, middle, and lower terciles); and (b) the top five AN and 
top five BN events. For the 40 years of our 1970-2009 study period, each tercile 
represents approximately 13 years (one-third of the 40 years). The resulting 
tercile categories for HOA PR were: AN PR is > 3.5 mm/day; NN PR is < 3.5 
mm/day but > 3.1 mm/day; and BN PR is < 3.1 mm/day. The resulting top five 
AN (BN) Oct-Nov HOA PR years were 1972, 1977, 1994, 1997, 2006 (1971, 
1974, 1998, 2005, 2007). 

We then developed composites of global and regional anomalous PR, 
SST, GPH, and other climate variables (Table 1) for the top five AN and top five 
BN events. We also composited according to the top five AN and top five BN 
events for the three climate variation indices (MEI, Nino3.4, and DMI) and for the 
potential predictors. We used these composites to assess the patterns and 
processes during extreme event years. In particular, we used the extreme event 
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composites to identify potential predictors and the dynamical processes that 
influence Oct-Nov precipitation in the predictand region. We focused in 
composites based on just the top five AN and BN events, because composites 
based on more than five events tend to include weaker events and tend to 
obscure important climate patterns, processes, and mechanisms of interest 
(LaJoie 2006). All the anomalies in our study were calculated using a base 
period of 1968-1996 (cf. ESRL 2010). 


Table 1. Climate system variables that we analyzed to identify: (a) potential 
predictors of HOA PR; and (b) the physical processes that lead to climate 
variations in HOR PR. All variables were assessed on a global scale, 
except for soil moisture. 


Variable 

Abbreviation 

Units 

Sea surface temperature 

SST 

degrees Celsius rc) 

Air temperature at 10 tn 

T air air temperature 

degrees Celsius ("C) 

200 hPa geo potential heights 

200 hPa GPH 

geopotential meters <GP Ml 

850 hPa geopotential heights 

850 hPa GPH 

geopotential meters <GP M| 

200 hP a vector Winds 

- 

meters second (m s) 

850 hP a vector winds 

- 

meters second (m s) 

Specific humidity 

SH 

grams centimeters (g cm3) 

500 hPa vertical velocity 

Omega 

Pascals second (Pas) 

P recipitation rate 

PR 

millimeters day (mm day) 

Outgoing longwave radiation 

OLR 

watts meter2 (w m2) 

Volumetric Soil Moisture (0-10cm) 

soil moisture 

fraction 


We calculated bi-monthly (e.g., Oct-Nov) correlations at lag times of zero 
to six months between the predictand PR time series and the climate system 
variables shown in Table 1. We also calculated similar correlations between the 
three climate variation indices and global PR. These correlations were analyzed 
to identify the spatial patterns and teleconnections associated with climate 
variations in HOA PR. For correlations based on the full study period of 1970- 
2009, correlations greater than +/- 0.304 were considered statistically significant 
at the 95% confidence interval based on the standard normal distribution of a two 
tailed Mest (Wilks 2006). When the dataset was reduced to 10 years and 15 
years for our optimum climate normal analyses, correlations greater than +/- 0.48 
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and +/- 0.57, respectively, were considered statistically significant at the 95% 
confidence interval. The correlations were conducted using Excel and the ESRL 
web sites (ESRL 2010). 

3. Predictor Selection 

In our study, we defined potential predictors as variables with significant 
long lead correlations and dynamically plausible teleconnections with the 
predictand. All the predictors we analyzed represented area averaged values 
(e.g., the average of SST in a 5 x 5 degree region). Multiple potential predictor 
variables were investigated (see Table 1), primarily via correlations with the 
predictand. SST at several locations was found to have the consistently highest 
significant correlations with the predictand at 0-6 month lead times. The method 
for selecting the predictor region was the same as that used for selecting the 
predictand region (see Chapter II, Section B.1). After selecting the predictors, we 
conducted composite analyses of the top five AN and top five BN predictor 
events to identify the patterns and dynamical processes associated with extreme 
predictor events and their relationships to extreme predictand events. 

4. Predictor and Predictand Time Series 

After choosing the predictors, we compared the predictor and predictand 
time series at all lead times. This allowed us to: (a) do visual quality control 
checks of the predictor-predictand correlations; (b) identify interannual to decadal 
variations in the relationships between the predictors and predictand; (c) identify 
periods and lead times for which the correlations were relatively strong and 
weak; and (d) identify case studies for additional composite, correlation, and 
dynamical analyses, and hindcast testing. 

5. Long-Range Forecasting Methods 

We developed and assessed several LRF methods for this study in order 
to improve upon: (a) the LTM forecast methods commonly used by DoD 
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forecasters, and (b) the LRF methods presently being used by 14 WS and other 
non-DoD organizations. A number of prior studies have shown that climate 
analyses and forecasts that account for deviations from LTMs have several 
potential benefits to DoD planners by providing more accurate and reliable 
assessments and forecasts than those based on LTMs alone (LaJoie 2006; 
Vorhees 2006; Hanson 2007; Moss 2007; Turek 2007; Crook 2009; Heidt 2009; 
Ramsaur 2009). In addition, LRFs based on customized predictors selected for a 
specific forecast target often have more skill than those based on general 
predictors, such as ENLN indices (e.g.. Crook 2009; Heidt 2009). 

The deterministic and probabilistic LRF methods we used in this study 
include tercile matching, averaged linear regression, non-averaged linear 
regression, composite analysis, and optimal climate normals. We applied these 
forecast methods to several predictor-predictands pairs and to several hindcast 
and forecast time periods. The main results from each of these methods were 
categorical forecasts; that is, forecasts of the predictand tercile category—AN, 
NN, or BN HOA PR. 

A primary motivation for using such categorical forecast targets is that 
they are simpler and larger forecast targets than discrete value forecast targets, 
which tends to: (a) make the development of the forecast method simpler; (b) 
increase predictability at long lead times (e.g., leads of several months or 
seasons); and (c) simplify forecast verification (cf. van den Dool 2007). The 
disadvantage of using a categorical predictand is that the forecast gives only a 
range of predictand values and does not specify where within that range the 
actual predictand variable will occur. However, this disadvantage can be 
mitigated by selecting smaller categories than terciles (e.g., quintiles). This 
disadvantage may not be as important to end users of the forecasts as they are 
to the forecasters. This tends to be the case for operational planners who only 
need to know the likelihood that an environmental variable will be above or below 
a given operational threshold. The need for such threshold forecasts is common 
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for DoD planners, especially when planning operations at long lead times. 
Finally, if LRFs of more specific forecast targets do not have skill, then 
categorical LRFs may be the only viable option. 

All of the forecast methods use just one to three predictors, and are thus 
unlikely to predict (or explain) all of the predictand variations. However, by 
careful selection of the predictors (and the predictand), these methods have the 
potential to predict many of the predictand variations, and thereby produce skillful 
LRFs (cf. Hanson 2007; Moss 2007; van den Dool 2007; Heidt 2009). 

To assess the viability of our LRF methods, we used the methods to 
conducted hindcasts for the period 1970-2009. We then assessed the skill of 
these multi-decadal hindcasts according to the predictor, lead time, and forecast 
method. 

For all the forecast methods, we used Excel to help analyze the predictor 
and predictand data, generate the forecasts, and assess the performance of the 
forecasts. The methods are described in the following subsections. 

a. Method 1 - Tercile Matching 

For the tercile matching method, we then used the sign of the 
correlation between the predictor and predictand to determine the expected 
relationships between the predictor tercile categories and the predictand tercile 
categories. For example, if the correlation was negative, then the expected 
relationships would be: If AN (BN, NN) predictor, then expect BN (AN, NN) 
predictand. The forecast is based on the expected relationship between the 
known tercile category for the predictor and the unknown tercile category for the 
predictand. So, for the negative correlation example just given, an AN (NN, BN) 
predictor value would lead to a forecast of a BN (NN, AN) predictand value. To 
apply this method, we grouped the predictor and predictand seasonal values into 
their respective AN, NN, and BN tercile categories. We then determined for each 
year whether the predictor-predictand pair categories matched the expected 
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relationship (e.g., AN predictor and BN predictand, for the previous example), 
were opposite to the expected relationship (e.g., AN predictor and AN predictand, 
for the previous example), or were otherwise (e.g., AN predictor and NN 
predictand, for the previous example). For the years in which the pairs matched 
the expected relationship, the forecast was deemed accurate. For example, if we 
found a negative correlation between an SST predictor in Apr-May and HOA PR 
in Oct-Nov, then we used the occurrence of an AN (NN, BN) SST in Apr-May to 
produce a forecast of BN (NN, AN) FlOA PR for the following Oct-Nov. 

The advantages of the tercile matching method include its conceptual 
and computational simplicity. In principle, it can be thought of as a very 
simple categorical version of linear regression. 

b. Method 2-A veraged Linear Regression 

We used two linear regression methods to model the relationships 
between HOA PR and several predictors, based on the R1 dataset. Linear 
regression (LR) refers to any approach to modeling the relationship between a 
predictor variable and two or more predictand variables, such that the model 
depends linearly on unknown parameters that are estimated from predictor and 
predictand datasets (cf. Wilks 2006). If the relationship is linear, LR can be a 
useful method for long-range forecasting giving an exact value. The LR model 
does not provide a determination that there is a causal relationship between the 
predictor and predictand. However, the LR model may help corroborate other 
evidence of a causal relationship. 

Liner regression result in the derivation of a regression equation 
that has the form (cf. Wilks 2006): 

Predicted value = slope * predictor + y-intercept 

This equation is the LR model and is determined by analysis of the 
predictor and predictand datasets. The resulting forecasts are discrete. 
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deterministic forecasts of the predictand that we then assigned to the AN, NN, or 
BN tercile categories, with the terciles being defined by the predictand time 
series derived from the R1 dataset. 

The 1970-2009 hindcasts for the two LR methods were formulated 
and tested by using a leave-one-out cross validation approach (Wilks 2006). 
That is, data for the year to be hindcasted were withheld when developing the 
regression model that was applied to that year, with the model for that year being 
based on data from all the other years (cf. Jury 1996; Wilks 2006). New AN, NN, 
and BN terciles were determined from the remaining years, which were used to 
determine the tercile of the new hindcasted value. For example, to formulate the 
1970 PR hindcast, we excluded the SST and PR values for 1970 from the R1 
dataset, leaving the 1971-2009 SST and PR values to form new AN, NN, BN 
terciles, and the slope and y-intercept. For the 1971 PR hindcast, we excluded 
the 1971 SST and PR from the R1 dataset, leaving 1970 and 1972-2009 to form 
the new terciles, slope and y-intercept, and so on through 2009. 

For the average linear regression method, all of the slope and y- 
intercepts were averaged together to come up with a single slope and y-intercept 
regression equation model in which the predictor value for the excluded year was 
inserted into the equation to get a hindcasted PR for the excluded year. The 
hindcasted value was then assigned to one of the new tercile bins and compared 
to the actual PR for that year. The average LR method yields one regression 
model for the entire dataset. The LR models that use non-averaged slope and y- 
intercept values were very similar from year to year, and very similar to the 
average LR model, as might be expected when only one out of a total of 40 years 
is being omitted per hindcast year. These similarities indicated that the average 
LR model was a simple but reasonable alternative to the non-averaged LR 
model. These similarities also indicated that the non-averaged LR models were 
not overfit (i.e., were not overly reliant on predictor-predictand relationships that 
were unique to a few years and unrepresentative of most years). This result 
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helped validate the use of non-averaged LR models. We developed averaged 
regression equations for lead times of zero to six months. 


c. Method 3 - Non-A veraged Linear Regression 

The development of the non-averaged LR models is described in 
Section 5.b, but without the averaging of each year’s slope and y-intercept values 
to produce a single average LR model. So, in Method 3 we did not average the 
40 regression equations into one. For example, to hindcast the 1970 PR, we 
used the leave-one-out cross validation approach described in Section 5b, 
excluding 1970 from the dataset and produced a slope and y-intercept from the 
remaining 1971-2009 data. This gave us our regression equation from that 
dataset and became our LR model to hindcast the 1970 PR, and so on through 
2009. The difference between the predictions from the averaged and non- 
averaged regression equations were small, although sometimes large enough, to 
change the predicted PR to a different tercile (e.g. BN to NN, or NN to AN). We 
developed regression equations for all hindcast years 1970-2009 and at lead 
times of zero to six months. The non-averaged LR model based on the years 
1970-2009 was our preferred LR model for producing LRFs for 2010. Likewise, 
a non-averaged LR model for the years 1970-2010 would be our preferred LR 
model for LRFs of 2011. 

d. Method 4 - Composite Anaiysis Forecast 

We used the composite analysis forecast (CAF; MetEd 2010) 

process to generate probabilistic, categorical, long-range forecasts based on the 

historical occurrence, or conditional probability, of a certain predictand category 

occurring given the occurrence of a specified predictor category (e.g., the 

probability of AN PR in Oct-Nov given the occurrence of BN SST conditions near 

New Zealand in Apr-May). In the CAF method, the statistical significance of the 

predictor-predictand relationship needs to be determined, with a 90 percent 

confidence level being the standard for application of the method. The statistical 

significance: (a) is a standard measure of the confidence in the rejection of the 
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null hypothesis that there is no relationship the predictor and predictand; and (b) 
provides insight into whether the use of conditional composite probabilities based 
on the predictor-predictand relationships is valid in the development of forecasts 
(MetEd 2010). If there are no statistically significant results for any given lead 
time and predictor condition, then the CAP method would not be used. 

A forecast based on the CAP method involves three equations to 
compute the probability of a predictand state (AN, NN, or BN). Pigure 12 shows 
the three equations used to forecast the probability of HOA PR in the AN, NN, 
and BN categories (the left hand side terms in the equations) based on the 
relationship between the PR predictand and a SST predictor. The first terms on 
the right hand sides of the equations represent the conditional probabilities. Por 
example, the first right hand term in the first equation represents the probability of 
AN PR given the occurrence of AN SST. The conditional probabilities are 
determined from observations (e.g., the R1 dataset in our case). The second 
terms on the right hand side are the probabilities of the SST categories occurring. 
Por example, the second right hand term in the first equation represents the 
probability of AN SST occurring. In the standard application of the CAP method 
(MetEd 2010), the predictor probabilities are forecasted (e.g., forecasts of 
Nino3.4). 



Pigure 12. CAP method equations for forecasting the probability of AN, NN, 
and BN PR (the predictand categories; the terms on the left hand 
sides of the equations) given: (a); the conditional probabilities of the 
PR predictand categories given the occurrence of the specified 
SST predictor categories (the first terms on the right hand sides of 
the equations); and (b) the probability of SST predictor categories 
occurring (the second terms on the right hand sides). See main 
text for more information. The forecast lead time is not specified in 
these equations. Adapted from MetEd (2010). 
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The main steps in the CAF methods are: 

(1) Select predictand region of interest (e.g., HOA). 

(2) Select analysis and forecast periods and lead times (e.g., 1970- 
2009 dataset for predictand and predictor, Oct-Nov forecast target period, 0-6 
month leads). 

(3) Select predictor (e.g., SST). 

(4) Determine predictand and predictor categories (e.g., terciles). 

(5) Determine conditional probabilities based on datasets (e.g., AN 
PR in Oct-Nov occurred 65% of the time when BN SST conditions occurred near 
New Zealand in May-Jun). 

(6) Conduct risk analysis using hypergeometric distributions to test 
for statistical significance in predictor-predictand relationship. 

(7) Calculate composite analysis forecast for statistically significant 
relationships at a 90% confidence level. Use alternate LRF method is 
relationships are not significant. 

(8) Verify CAF method with hindcasts. 

The strength of the CAF method is that it exploits the same 
historical observed data used to generate LTMs, yet does so in a way that yields 
potentially enhanced forecast and decision-making information. The method 
does however depend on the availability, accuracy and representativeness of the 
predictor and predictand data. See MetEd (010) for more details on the CAF 
method. See Flansen (2007), Moss (2007), Crook (2009), and Fleidt (2009) for 
reviews and applications of the CAF method. 

In our study (as in Crook 2009 and Fleidt 2009), we chose to make 

the predictors non-forecasted variables (rather than forecasted ones, such as 

Nino3.4 or DMI). We did this because: (a) our predictors were methodically 

selected to maximize the predictor-predictand correlations at long lead times; (b) 

forecasts of our selected were not available; and (c) the selection of leading 
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predictors meant that once the predictor was observationally determined at the 
appropriate lead time, there was no need to forecast the predictor. Thus, in our 
modification of the standard CAF equations (Figure 12), the predictor probability 
is one for the predictor in the observed tercile and zero for all the other predictor 
probabilities (Figure 13). The modified CAF equations in Figure 13 represent the 
equations for a case in which the SST predictor at the forecast lead time is 
observed to be AN. Thus, for a CAF at this lead time, the SST predictor 
probability is one for AN SST and zero for NN and BN SST. This means that 
many of the terms in the standard CAF equations go to zero in the modified CAF 
equations. 


^jjAN[PR] 

1 AA?[55r] 

^ pBN[PR] 
1 AA^[55r] 


Figure 13. Equations used in our modified CAF process to calculate long- 
range forecasts of PR based on observed AN SST predictor 
conditions at the forecast lead time. The lead time is unspecified in 
these equations. Compare to the standard equations in Figure 12. 





A key advantage of our adaptation of the CAF method is that CAFs 
can be generated for a wide range of customized predictors, rather than just the 
few predictors that are forecasted at long leads (e.g., Nino3.4). This allows the 
use of predictors that are well correlated with the predictand but not forecasted. 
However, for operational use, the predictors must be well observed in real time, 
so that observed values are available at the appropriate lead times. Thus, for 
example, well-correlated SST predictors may be good predictors for the adapted 
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CAF method because SST observations are routinely available in real time. But 
well correlated near-surface ocean temperatures may not be good predictors, 
because near-surface ocean temperature observations are generally not 
available in real time (except for a few special locations (e.g., the equatorial 
Pacific and other special locations with moored buoy observations of ocean 
temperature). 


e. Optimal Climate Normals (OCN) Approach 

LTM climate conditions are the average of conditions observed 
over a number of years (e.g., 30 years; WMO 2010b). In a stationary climate, 
LTMs are good estimates of the upcoming expected value for a given variable 
(e.g., precipitation rate). However, if the climate system varies over periods of 
several years (e.g., 5-15 years), then taking the average of the most recent 
years may provide a better estimate of the expected value than would be 
obtained from a longer period mean (e.g., a better estimate than that obtained 
from the WMO standard 30 year mean) (van den Dool 2007). The standard way 
of applying the optimal climate normals (OCN) approach is to calculate a mean 
value for the predictand from a selected set of consecutive recent years, and 
then persist that mean to produce a forecast for the year that immediately follows 
the most recent year of the set (van den Dool 2007). CPC uses OCN for many of 
its seasonal to interannual forecasts, such as U.S. surface air temperature, 
precipitation, and North Atlantic hurricane activity (Livezey 2006; van den Dool 
2007). CPC uses the average of the most recent 10-15 years when generating 
operational OCN forecasts. 

For this study, we applied the OCN concept by using two 
approaches to account for recent relationships between the predictor and 
predictand. In the first, or basic, OCN approach, we used the same predictor that 
we used when hindcasting all 40 years of our 1970-2009 study period, but we: 
(a) determined our tercile categories for the predictand and predictor using just 
the most recent 10- and 15-year periods (2000-2009 and 1995-2009); (b) used 
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the 10- and 15-year tercile categories in generating hindcasts for those recent 
years using our four forecast methods; and (c) used the 10- and 15-year tercile 
categories in verifying the hindcasts. We termed this the basic OCN approach. 

In the second, or advanced, OCN approach, we conducted all our 
analyses (e.g., selection of predictors, composites, correlations, teleconnections) 
using data from just the most recent 10- or 15-year periods during 1970-2009 
(e.g., from just 1995-2009). From this, we identified two new, well correlated, 
and dynamically plausible predictors that we then used in the four forecast 
methods described in Sections 5.a-d, above. The resulting hindcasts were then 
verified against the observed conditions during the most recent 10 to 15 years. 
This allowed us to determine whether forecast methods that were based on 
conditions during the most recent 10-15 years could lead to improved forecast 
skill. We termed this the advanced OCN approach. 

In neither OCN approach did we replicate the standard OCN 
process in which the forecast is based on persisting the mean of the recent 10- 
15 years (van den Dool 2007). But, both of our approaches apply the basic OCN 
concept of emphasizing recent conditions when developing forecast methods 
and generating forecasts. We applied the OCN approach to all four of our 
forecast methods. Thus, we refer in this report to OCN approaches, rather than 
OCN methods. We distinguished our two OCN approaches and our other 
approach involving all of our study period years (1970-2009) by referring to them 
as the basic OCN approach, the advanced OCN approach, and the all years 
approach, respectively 

6. Hindcast and Forecast Verification Methods 

All of the predicted PR values from the different forecast methods, with 
and without applying OCN methods, were binned into AN, NN, and BN terciles to 
generate categorical PR hindcasts for 1970-2009. These tercile category 
hindcasts were then compared to the tercile category of the actual PR to assess 
the skill of the forecast methods. This verification was done using a 3x3 
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contingency table method to determine the number of hits, misses, false alarms 
(FA), and correct rejections (CR) (Wilks 2006). The 3x3 table accounts for the 
three predictor terciles and the three predictand terciles. This contingency table 
method was applied separately for each predictor, each forecast category (AN, 
NN, BN), each forecast method, and each lead time. Table 2 shows the 
assignments that would be made to a 3x3 contingency table for a forecast of AN 
PR based on BN SST predictor conditions. The four types of assignments are; 

1. Hit: event forecasted and observed 

2. Miss: event not forecasted but observed 

3. Faise aiarm (FA): event forecasted and was nof observed 

4. Correct rejection (CR): event nof forecasted and was not observed 

For the case shown in Table 2: a hit would be assigned when AN PR was 
forecasted and observed; a miss would be assigned when AN PR was observed 
but not forecasted; FA would be assigned when AN PR was forecasted but not 
observed; and CR would be assigned when AN PR was not forecasted and was 
not observed. 

Table 2. Example of contingency table assignments for the verification of 
forecasts of AN PR when using a categorical SST predictor, for the case 
in which AN SST leads to a forecast of BN PR. 



Predictand 

Predictor 

AN PR 

NN PR 

BN PR 

AN SST 

Miss 

Correct Rejection 

Correct Rejection 

NN SST 

Miss 

Correct Rejection 

Correct Rejection 

BN SST 

Hit 

False Alarm 

False Alarm 


Table 3 shows the hit and miss assignments for AN, NN, and BN PR for a 
case in which there is a negative correlation between the SST predictor and the 
PR predictand. 
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Table 3. Schematic contingency table of hit and miss assignments based on 
a negative correlation between a SST predictor and a PR predictand. 



Predictand 

Predictor 

AN PR 

NN PR 

BN PR 

AN SST 

Miss 

Miss 

Hit 

NN SST 

Miss 

Hit 

Miss 

BN SST 

Hit 

Miss 

Miss 


With the exception of the tercile matching method, each LRF method 
yields a discrete (non-categorical) forecast of PR or of the probability of a PR 
category. In order to apply the contingency table method to these discrete and 
probabilistic forecasts, we: (a) assigned the discrete forecasts to the terciles 
categories based on the PR terciles determined from the R1 dataset; and (b) in 
some cases, assigned the probabilistic forecasts to tercile categories based on 
the highest probability category. This allowed the forecasts from all the forecast 
methods to be readily compared. Examples of predictand and predictor terciles 
are shown in Figure 14. 
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HOA PR Predictand 



Figure 14. Time series of Oct-Nov HOA PR predictand (top panel) and New 
Zealand SST predictor (bottom panel) for 1970-2009. of the bold 
black and red horizontal lines denote the boundaries between the 
AN, NN, and BN categories for the predictand and predictor, 
respectively. For example, PR values are categorized as AN when 
> 3.5 mm/day; NN when < 3.5 mm/day but > 3.1 mm/day, and BN 
when < than 3.1 mm/day. 
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The contingency table assignments are used to calculate several forecast 
verification metrics, including: percent correct (PC) or accuracy rate, FA ratio 
(FAR), probability of detection (POD), and Heidke skill score (MSS) (Wilks 2006). 
These metrics were calculated for each predictor, each forecast category (AN, 
NN, BN), each forecast method, and each lead time. In addition, to see the 
variation in skill scores among the lead times and predictor events, we 
calculated: (a) category-average skill scores averaged over all categories (AN, 
NN, and BN) for each lead time; and (b) total skill scores averaged over all 
categories and all lead times. 

The following equations show how each metric was calculated (Wilks 

2006): 

1. PC: = (hits + CRs) / total 

2. FAR: = FAs / (hits + FAs) 

3. POD: = hits / (hits + misses) 

4. HSS: = (hits + CR - expected correct) / (total - expected correct) 

Here, the term total refers to the total number of forecasts. The HSS 
helps answer the question: “What was the accuracy of the forecast relative to 
that of random chance?” where random chance is determined from LTM 
climatology. HSS is based on a scale from one to minus one—one being the 
best score and minus one being the worst. Having a HSS of less than zero 
(zero) means that using the chosen predictor-predictand pair as your forecast 
method is worse (no better) than using LTM climatology. 

The following three metrics thresholds were used to identify forecasts with 
relatively high skill, especially to determine the viability of each predictor- 
predictand pair (e.g., their potential viability for use in the CAF method). 

(1) PC greater than 50% 

(2) POD equal to or greater than FAR 

(3) HSS values greater than 0.3 
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Forecasts and pairs for which all three thresholds were met were identified 
as especially skillful and viable, respectively. Forecasts and pairs that did not 
meet these thresholds were rejected from further consideration. 

C. SUMMARY OF FORECAST METHODS 

Figure 15 summarizes the major data sets and methods we used to 
produce long-range hindcasts and forecasts. 



Choose Predictor 


Develop Method 4: Composite 
Analysis Forecast 

1. Perform composite analysis 

2. Condjct risk analysis and 
determine statistical significance 
Are CAP results statistically 
significant (> 90 %) 


i*oo' 


I 



Long Range Method Development & Testing 

1. Deterministicmethods 

a) Method 1: tercile matching 

b) Method 2; avgd linear regression 

c) Method 3: non-avgd linear regression 

2. Hindcasts, verification & forecasts 

Is predlctor/predlctand pair viable? 


Base LRF on 
deterministic or 
probabilistic forecast 


73Er 

I I 


Evaluate other 
potential 

predictor variables 


Base LRF on LTM and/or conditional means 


1 

Provide tailored 


long lead forecast 

w 

support products 


Figure 15. Flow chart showing main datasets and methods used to conduct 
climate analyses and long-range forecasts in this study. 
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III. RESULTS 


A. CLIMATE ANALYSES 


1. Composites Conditioned on Extreme PR in the HOA 


We selected our predictand as Oct-Nov HOA PR in the region bounded 
by 5S-10N, 30-55E. This predictand region is shown by the black box in Figure 
16 and includes almost all of the broad HOA region (see Chapter I, Section G.4), 
with the exception of northern Ethiopia, Eritrea, Djibouti, and most of Sudan. We 
chose this as our predictand region because of the tendency for similar spatial 
and temporal patterns of variability in Oct-Nov PR for the points within this 
region, based on 1970-2009 PR from the R1 dataset (see Chapter II, Section 
B.1). 



Oct to Nov: 



Figure 16. The Oct-Nov HOA PR predictand region for our study (black box) 
overlaid on the LTM Oct-Nov PR (color shading; mm/day; cf. 
Figure 10). PR image created using ESRL (2010). 
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The time series of Oct-Nov PR for 1970-2009 for the predictand region is 
shown in Figure 17. Note the large interannual variations, especially those 
exceeding plus or minus one standard deviation. This PR times series was used 
to: (a) identify the top AN and BN events during 1970-2009 for use in 
constructing composite anomalies; and (b) to conduct correlation and 
teleconnection analyses (see Chapter II, Section B.2). Figure 17 also shows the 
MEI and DMI times series, which reveal a general correspondence between; (a) 
ENLN and the lOD (Twigg 2007); and (b) the predictand PR and ENLN and the 
lOD (e.g., Ogallo 1988; Black at al. 2002). Note in particular the general 
tendency for AN (BN) FlOA PR in Oct-Nov during EN (LN) and positive 
(negative) lOD. These general correspondences suggest that: (a) Oct-Nov HOA 
PR is teleconnected to global scale climate variations; and (b) ENLN and lOD 
indices (e.g., the MEI and DMI) may be useful predictors of Oct-Nov HOA PR. 
However, Figure 17 shows some periods in which one or more of these general 
correspondences did not hold (e.g., 1980 and a number of years during 1984- 
1993). In addition. Figure 15 does not show how the ENLN and lOD indices are 
related to the HOA PR when the indices lead the HOA PR by one or more 
months. So, it is unclear from Figure 17 alone just how useful these indices 
would be as long-range predictors of Oct-Nov HOA PR. 
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Figure 17. Time series of HOA PR (predictand, purple line, mm/day), and El 
Nino-La Nina (ENLN) and Indian Ocean Dipole (lOD) indices (MEI, 
red line, and DM I, green line, respectively) for Oct-Nov 1970-2009. 
Note the similarities in the interannual variability of the predictand 
and the two indices. Indications of teleconnections between HOA 
PR, ENLN, and lOD are highlighted during: (a) a strong LN and 
negative lOD period in Oct-Nov 1974-75; and (b) a strong EN and 
positive lOD period in Oct-Nov 1997. Relationships such as these 
suggest the usefulness of evaluating ENLN and lOD indicators as 
potential predictors for HOA PR. 

To further investigate the regional and global scale patterns and 
processes associated with HOA PR variations, we constructed conditional 
composites of the anomalies during the Oct-Nov periods in which our HOA 
predictand region had the highest and lowest PRs (i.e., the most extreme AN and 
BN periods). We designated these as the wet and dry composites, respectively. 
See Chapter II, Section B.2 for more on what years these were and how they 
were identified. These and similar composites conditioned on the occurrence of 
extreme AN and BN PR in our predictand region helped us identify regional and 
global scale patterns and processes (e.g., teleconnection) that are related to PR 
extremes in our predictand region. 

Figure 18 shows the regional and global wet and dry composites for PR. 
Note that AN (BN) PR in the HOA is associated with: (a) AN (BN) PR in a large 
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region that includes the HOA and extends over the western 10; (b) BN (AN) PR 
over tropical Africa to the west of the HOA; (c) BN (AN) PR over the eastern 
tropical 10, maritime continent (MC), and South Pacific Convergence Zone 
(SPCZ); and (d) AN (BN) PR over the central and east tropical Pacific. These 
patterns are consistent with the results shown in Figure 17 and with prior studies 
of teleconnections that affect east African and HOA PR (e.g., Ogallo 1988; 
LaJoie 2006; Vorhees 2006). They are also very useful in identifying: (a) the 
dynamical processes that cause climate variations in HOA PR; (b) potential 
predictand regions; and (c) potential predictors of HOA PR. 


Regional View Global View 



-3 -2.6 -2-2 -1.8 -1.4 -1 -0.6 -0.2 0.2 0.6 ! 1.4 1.8 2.2 2.6 3 


Figure 18. Conditional composite anomalies of Oct-Nov precipitation rate (PR; 

mm/day) for the five most extreme AN (wet) and BN (dry) events of 
HOA PR during Oct-Nov for 1970-2009. The left (right) panels 
show regional (global) scale views. See Chapter II, Section B.2 for 
the five most extreme AN and BN PR event years. The red box 
outlines our HOA predictand region. Note the generally opposite 
patterns in the wet and dry composites, and the correspondence to 
PR patterns associated with ENLN and the lOD. PR anomaly 
patterns such as these suggest possible predictors and linkages 
associated with the climate variations in the HOA. Images created 
using ESRL (2010). 
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Figure 19 shows the global scale wet and dry composite SSTAs, with AN 
(BN) PR in the HOA being associated with SSTAs that are: (a) positive (negative) 
in the western lO; (b) negative (positive) in the eastern tropical lO, western 
tropical Pacific, SPCZ-New Zealand region, and in and northeast of the 
Philippine Sea; and (c) positive (negative) in the central-eastern tropical Pacific 
and along the west coast of the Americas. These patterns are consistent with 
the PR anomalies shown in Figure 19, and with the relationships between FlOA 
PR and ENLN and lOD indicated by Figures 18 and 17. The approximately 
opposite patterns in the wet and dry composites suggest that the global scale 
dynamical processes associated with interannual variations in HOA PR are 
predominantly linear. 



Figure 19. Conditional composite anomalies of sea surface temperature (°C) 
for five most extreme wet (left panel) and dry (right panel) HOA PR 
events during Oct-Nov 1970-2009. Note in the wet (dry) 
composite the positive (negative) SSTAs in the WIO and central- 
eastern tropical Pacific and the negative (positive) SSTAs in the 
tropical ElO and western tropical Pacific, consistent with positive 
(negative) lOD and EN (LN) conditions. Images created using 
ESRL (2010). 


Figure 20 shows the wet and dry composite anomalies for the 200 hPa 
geopotential height (GPHs). These anomalies are useful in determining the 
large-scale atmospheric circulation anomalies associated with wet and dry HOA 
PR periods. The anomalies over the tropics and subtropics show indications of 
tropical Rossby-Kelvin wave responses to the SST and precipitation anomalies 
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shown in Figures 18 and 19 (e.g., the paired anticyclones over Africa and the 
central-eastern tropical Pacific in the wet composite; the negative (positive) GPH 
anomalies over subtropical Asia and subtropical southern 10 in the wet (dry) 
composite). They are also consistent with prior studies of the Rossby-Kelvin 
wave response to ENLN and lOD (e.g., LaJoie 2006; Vorhees 2006). 



Figure 20. Conditional composite anomalies of 200 hPa GPH (gpm) for five 
most extreme wet (left panel) and dry (right panel) HOA PR events 
during Oct-Nov 1970-2009. Note in the wet (dry) composite the 
positive (negative) GPHAs over Africa and the central-eastern 
tropical Pacific (central and southern Africa and central-eastern 
tropical Pacific), consistent with positive (negative) lOD and EN 
(LN) conditions. These features are indicative of Rossby-Kelvin 
wave responses to anomalous convection over the WIO, MC and 
central-eastern tropical Pacific. Images created using ESRL 
( 2010 ). 


In the wet composite, above normal 200 hPa GPHs exist over much of 
Africa and eastward along the equator in response to enhanced convection over 
east Africa and the western 10 (cf. Figure 18; LaJoie 2006; Vorhees 2006). At 
about 30 degrees to the north and south of the MC, below normal GPH 
anomalies occur in response to suppressed convection over the MC. The dry 
composite GPH anomalies are nearly opposite. This is most visible in the above 
normal GPH anomalies over southern Asia and Australia extending westward 
over northern Africa and southern 10, respectively, as a response to the 
enhanced convection in the MC. In both the wet and dry composites, the 200 
hPa GPH anomalies over the equatorial Pacific are consistent with the ENLN 

related PR and SST anomalies in this region (Figures 18-19). The 200 hPa 

54 




































GPH anomalies are more pronounced in the wet composite, possibly because 
the HOA PR, ENLN, and lOD variations in the wet years were more extreme 
deviations from their LTMs than those in the dry years (cf. Figure 16). The wet 
and composite shows clear indications of an anomalous wave train extending 
along an arching path from the maritime continent (MC) sector into the southern 
hemisphere extratropics, reaching its most poleward latitude at about 90° W, and 
then extending back into the tropics over the South Atlantic. There is less clear 
evidence in the dry composite of an oppositely signed anomalous wave train in 
the southern hemisphere. 

The corresponding GPH composite anomalies at 850 hPa (not shown) 
are, in the tropics and subtropics, nearly opposite to those at 200 HPa (Figure 
20). In particular, in the wet (dry) composite, the 850 hPa GPH anomalies are 
below (above) normal over eastern Africa and the WIO, and above (below) 
normal over the MC. These GPH anomaly results are consistent with those at 
200 hPa and with those associated with ENLN and the lOD (Rassmusson and 
Carpenter 1982; Saji et al. 1999; LaJoie 2006; Vorhees 2006). 

The 850 hPa zonal wind anomalies in the wet composite (Figure 21, left 
panel) clearly shows anomalous easterlies from the MC converging with 
anomalous westerlies across central Africa in and near the HOA. The 
anomalously convergent low level wind flow into the HOA aids precipitation 
development in the HOA by transporting energy and moisture into the region 
(and similarly for the central-eastern equatorial Pacific, and oppositely for the 
MC). The dry composite (Figure 21, right panel) shows opposite low level wind 
anomaly patterns, with low level convergence centered over the MC and 
divergence over the HOA and central equatorial Pacific. In addition, in the dry 
composite, northerly flow from the Caspian Sea region toward the HOA occurs, 
indicating a tendency for anomalously stable and dry conditions over the HOA. 
The wet (dry) composite shows anomalous low-level divergence (convergence) 
in the SPCZ region. 
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Figure 21. Conditional composite anomalies of 850 hPa zonal wind (m/s) for 
the five most extreme wet (left panel) and dry (right panel) HOA PR 
events during Oct-Nov 1970-2009. The red arrows represent the 
sense of the wind flow anomalies. Note in the wet (dry) composite 
the anomalous low-level convergence (divergence) over the HOA 
and divergence (convergence) over the MC, consistent with above 
normal PR in the HOA and with positive (negative) lOD and EN 
(LN) conditions. Images created using ESRL (2010). 


Figure 22 shows wet and dry composite anomalies for outgoing longwave 
radiation (OLR). The OLR anomaly (OLRA) patterns are very similar to those for 
PR and SST (Figures 18-19). This is not surprising, since small changes in 
SSTs can create large anomalies in the energy and moisture available for 
atmospheric convection. OLR is a good proxy for convection in the tropics, with 
positive (negative) OLRAs corresponding to negative (positive) convection 
anomalies. In the wet composite, the OLRAs are below (above) normal over 
eastern Africa and western 10 and extending northward into southwest Asia, 
including Afghanistan, and in the central-eastern equatorial Pacific (West Africa, 
MC, and SPCZ region). The opposite OLRAs occur in the dry composite. We 
point out the OLR anomalies in southwest Asia because of the ongoing military 
operations in Afghanistan. This topic has been investigated by Vorhees (2006), 
but additional research on long-range forecasting of precipitation in southwest 
Asia using the methods applied in our study should be considered for improving 
precipitation forecasts for operations in this region. 
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Figure 22. Conditional composite anomalies of OLR (W/m^) for five most 
extreme wet (left panel) and dry (right panel) HOA PR events 
during Oct-Nov 1970-2009. OLR is a good proxy for convection 
and precipitation in the tropics, with positive (negative) OLRAs 
indicating negative (positive) convection anomalies. Note in the 
wet (dry) composite the negative (positive) OLRAs in the HOA, 
WIO, and central-eastern tropical Pacific, and the positive 
(negative) OLRAs in the tropical eastern lO and western tropical 
Pacific, consistent with wet (dry) conditions in the HOA, and 
positive (negative) lOD and EN (LN) conditions. Images created 
using ESRL (2010). 


2. Correlations 

We correlated our Oct-Nov HOA PR predictand time series (Figure 17) 
with a number of other potential predictors for the period 1970-2009, including: 
(a) the MEI, DMI, other climate variation indices (e.g., indices of the North 
Atlantic Oscillation and Arctic Oscillation); and (b) global SST and the other 
climate variables listed in Table 1. We found that SST and SST based climate 
variation indices (e.g., MEI, Nino3.4, DMI) in certain regions were: (a) well 
correlated with the predictand, with the predictor leading the predictand by zero 
to eight months; and (b) had higher correlations than non-SST related variables 
at the longer lead times (e.g., four to six months). The correlations with SST in 
the subtropical southwest Pacific north of New Zealand were especially strong 
and were statistically significant out to eight months lead. However, we chose to 
investigate leads only out to six months to focus on the leads at which our LRFs 
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had the greatest skill. For our optimal climate normals investigations, we 
conducted additional correlation analyses for the most recent years in the 1970- 
2009 study period (see Chapter II, Section B.5.e, and Chapter III, Section C). 

In the correlation maps shown in this section (e.g.. Figure 23), spatial 
correlation patterns are shown for correlations in which our HOA PR predictand 
lags the variable of interest (e.g., SST) by zero to six months. Areas of positive 
(negative) correlations indicate that HOA PR tends to increase (decrease) zero to 
six months after the variable of interest increases (decreases). The significance 
levels for the different correlations conducted in our study are described in 
Chapter II, Section B.2. We generated correlation maps for all the variables 
listed in Table 1. However, we only show in this report the maps for correlations 
involving SST and HOA PR, since we determined that SSTs were more viable 
predictors than the other variables we considered. This was because overall and 
with respect to the other variables: (a) the correlations with SST were relatively 
strong, significant, and persistent at all leads, especially at longer leads; (b) the 
dynamical links between SST and our predictand were relatively straightforward 
at all leads; and (c) SST is relatively well observed in real time and thus could be 
a feasible predictor for operational LRFs. 

Figure 23 shows the correlations of global SST with the Oct-Nov HOA PR 
predictand at lead times of one to four months. Note the similarity in the 
correlation patterns to the wet and dry composite SSTA patterns (Figure 19). 
This indicates that: (a) the composite SSTAs are representative of the general 
temporal and spatial patterns associated with wet and dry events; and (b) the 
correlations provide a good representation of the major dynamical processes 
associated with extreme wet and dry events. Figure 23 shows an area centered 
near 30S, 180W (red box north of New Zealand) that has a high negative 
correlation with the predictand (blue box in panel a), with correlations > -0.5 at 
leads of zero to four months. These correlations indicate that SST in this region 
tends to decrease (increase) before and when HOA PR increases (decreases). 
These correlations, significant at the 95% confidence level, may indicate: (a) 
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potential mechanisms that drive climate variations in HOA PR; and (b) potential 
predictors for producing skillful long-range forecasts of HOA PR. 



Figure 23. Correlations of HOA PR in Oct-Nov with global SST for 1970-2009 
with SST leading by: (a) one month (Sep-Oct), (b) two months 
(Aug-Sep), (c) three months (Jul-Aug), and (d) four months (Jun- 
Jul). The blue box represents the Oct-Nov HOA PR predictand 
region. The red box indicates a region north of New Zealand of 
persistently strong and significant negative correlation between 
SST and the predictand. Note the similarities between the 
correlation patterns in this figure and the SSTA patterns in Figure 
18. Images created using ESRL (2010). 


Figure 23 shows other areas of significant correlations (e.g., areas of 
negative correlation centered near 20N, 165E and in the eastern tropical 10; an 
area of positive correlation in the central-eastern tropical Pacific). But, the area 
north of NZ has the largest correlations over the largest area at the longest lead 
times (out to eight months; not shown). 
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3. 


Predictor Selection 


Based on the correlation results, we investigated further the use of SST 
north of New Zealand (NZ) as a predictor of Oct-Nov HOA PR. To do so, we 
began by determining the specific region north of NZ to use as our predictor 
region, following the same basic process we used to determine our predictand 
region (see Chapter II, Section B.1). In this process, we evaluated five potential 
predictor regions within the NZ high correlation region and chose as our NZ 
predictor region the area bounded by 18-36S, 165E-165W. We chose this region 
over the other potential NZ regions because the SSTs in this region had: (a) 
significant correlations at lead times of zero to six months; and (b) the highest 
mean correlation when averaged over all months of all the potential NZ and 
global scale climate variation index predictors. We named the SST in this region 
the NZ SST predictor. 

Figure 24 shows the magnitude of the correlations between the predictand 
and NZ SST, MEI, Nino3.4, and DMI at leads of one to six months. Note that the 
NZ SST predictor: (a) is significantly correlated at all leads; (b) is better 
correlated to the predictand than the indices at all leads, especially at longer 
leads; and (c) is not as well correlated as the indices at a lead of zero months. 

Figure 24 shows that the ENLN (MEI and Nino3.4) and lOD (DMI) 
correlations dropped below the 95% significant level after four and five month 
lead times, respectively. The lower correlation for Nino3.4 at longer lead times 
may provide clues as to why the 14 WS LRFs for their HOA region have low skill. 
Of course, the results in Figure 24 are for our HOA region, not the 14 WS HOA 
region that extends further to the north than our region, and for a particular 
forecast valid period (Oct-Nov). Still, these results do suggest that while 
established climate variation indices that are well correlated with a predictand 
may be good predictors, they might not be the best possible predictors. 
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These results indicate the potential value in: (a) doing a methodical search 
for potential predictors, rather than choosing a pre-established climate variation 
index as a predictor; and (b) selecting a predictor based on results from multiple 
lead times, rather than selecting a predictor based mainly on zero lag 
correlations. 


Comparison of Potential Predictors 

0.7 



Season 

-•-NZ16SW -#-«£! - Nino34 -OMI 


Figure 24. Magnitude of correlations between Oct-Nov HOA PR predictand 
and four potential predictors of this predictand, with the predictors 
leading by zero to six months. The four predictors are; NZ SST, 
MEI, Nino3.4, and DMI. The colored text boxes on the right of the 
figure describe the climate variations associated with the last three 
of these potential predictors. Correlations > 0.304 indicate 
significance at the 95% level (dashed line). 


To further test the NZ SST predictor, we visually compared the 1970-2009 
predictand time series to the NZ SST predictor time series for all lead times. 
Figure 25 below shows these comparisons for leads of 0, 3, and 6 months. In 
general, the three predictor time series are in phase with each other (positively 
correlated), but with some notable exceptions (e.g., 2006). In general, the three 
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predictor time series are out of phase with the predictand time series at all lead 
times (negatively correlated). However, for 1983-1993 the predictor and 
predictand tended to be in phase (positively correlated). Some prior studies of 
correlations between lO SST with East African rainfall have found a similar 
decadal scale correlation pattern reversal (e.g., Clark et al. 2002). The reasons 
for such reversals need further investigation. The predictor-predictor and 
predictor-predictand comparisons in Figure 25: (a) help clarify why the predictor- 
predictand correlation is strong t multiple leads; and (b) indicate the periods for 
which hindcasts based on the NZ SST predictor may do well and poorly. 



Figure 25. Time series of HOA Oct-Nov PR (mm/day) predictand in Oct-Nov 
(black line) and NZ SST (°C) predictor in (a) Oct-Nov (red line), (b) 
Jul-Aug (purple line), and (c) Apr-May (blue line) during 1970- 
2009. Note the persistence in the SSTs indicated by the similarities 
in the SST time series. Note also the indications of negative 
correlations between the predictand and NZ SSTs. 


Figure 26 shows the regions for the predictor-predictand pair that we used 
in our initial LRF model development and hindcasting. 
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Figure 26. Correlations of HOA PR in Oct-Nov with global SST in Sep-Oct 
based on data from 1970-2009.. These and related correlations 
(Figures 23-24) led us to select SST north of New Zealand (red 
box) as the initial predictor for use in forecasting HOA PR in Oct- 
Nov (our predictand; blue box) at leads of one to six months. 
Image created using ESRL (2010). 


B. LONG-RANGE HINDCASTS: ALL YEARS APPROACH 

1. Tercile Matching Method with NZ SST Predictor 

We used the NZ SST predictor to generate hindcasts for 1970-2009 at 
leads of one to six months using all four forecasting methods (see Chapter II, 
Section B.5). This resulted in a large number of hindcasts. Representative 
examples of the results from these hindcasts are presented in this report, along 
with summaries of the results for all of the hindcasts. These summaries are 
based mainly on the verification methods described in Chapter II, Section B.6. 

Table 4 summarizes the results of the hindcasts made using forecast 
Method 1 (tercile matching) based on the NZ SST predictor. The overall hindcast 
performance was very good, especially for the AN and BN categories. For 
example, HSS for the AN and BN categories is consistently above zero, 
especially at shorter leads. The lower performance for the NN category is typical 
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of many LRF systems, due in part to the systems being mainly designed, as ours 
are, to forecast deviations from the NN category (i.e., deviations from LTM 
conditions). Overall, our performance criteria (Chapter II, Section B.6) were met 
(i.e., PC greater than .50, POD greater than FAR, HSS of 0.3 or greater). These 
results, and similar results from the use of the other forecasting methods, 
indicate that our NZ SST predictor is: (a) viable when used with forecasting 
Method 1; and (b) likely to be viable when used in our other forecasting methods. 
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Table 4. Verification results from hindcasts generated using Method 1 
(tercile matching). The hindcasts predicted Oct-Nov HOA PR during 
1970-2009 based on the NZ SST predictor, with the predictor leading by 
zero to six months. Columns A, B, C, and D represent the number of hits, 
false alarms, misses, and correct rejections, respectively, for predictions 
of AN, BN, and NN PR. Verification metrics consist of percent correct (% 
Corr), false alarm rate (FA Rate), probability of detection (POD), and 
Heidke skill score (HSS). See Chapter II, Section B.6 for more 
information on our verification methods. 



A 

B 

C 

D 

Verification Metrics 


Hits 

FA 

Misses 

Corr. Rej. 

% Corr 

FA Rate 

POD 

HSS 

Results for hindcasts at lead times of zero months 

AN PR 

8 

5 

6 

21 

0.7 

0.4 

0.6 

0.4 

BN PR 

9 

6 

4 

21 

0.8 

0.4 

0.6 

0.5 

NN PR 

6 

6 

7 

21 

0.7 

0.5 

0.5 

0.2 

Results 

or hindcasts at lead time of one month 

AN PR 

10 

5 

4 

21 

0.8 

0.3 

0.7 

0.5 

BN PR 

6 

4 

7 

23 

0.7 

0.4 

0.6 

0.3 

NN PR 

6 

9 

7 

18 

0.6 

0.6 

0.4 

0.1 

Results for hindcasts at lead time of two months 

AN PR 

12 

12 

2 

14 

0.7 

0.5 

0.5 

0.3 

BN PR 

1 

1 

12 

26 

0.7 

0.5 

0.5 

0.1 

NN PR 

5 

9 

8 

18 

0.6 

0.6 

0.4 

0.1 

Results for hindcasts at lead time of three months 

AN PR 

7 

7 

7 

19 

0.7 

0.5 

0.5 

0.2 

BN PR 

9 

7 

4 

20 

0.7 

0.4 

0.6 

0.4 

NN PR 

1 

9 

12 

18 

0.5 

0.9 

0.1 

-0.3 

Results forhindcasts at lead time of four months 

AN PR 

6 

5 

8 

21 

0.7 

0.5 

0.5 

0.2 

BN PR 

6 

6 

7 

21 

0.7 

0.5 

0.5 

0.2 

NN PR 

4 

13 

9 

14 

0.5 

0.8 

0.2 

-0.2 

Results for hindcasts at lead time of five months 

AN PR 

8 

6 

6 

20 

0.7 

0.4 

0.6 

0.3 

BN PR 

6 

7 

7 

20 

0.7 

0.5 

0.5 

0.2 

NN PR 

3 

10 

10 

17 

0.5 

0.8 

0.2 

-0.1 

Results ■ 

or hindcasts at lead time of six months 

AN PR 

8 

7 

6 

19 

0.7 

0.5 

0.5 

0.3 

BN PR 

6 

6 

7 

21 

0.7 

0.5 

0.5 

0.2 

NN PR 

3 

10 

10 

17 

0.5 

0.8 

0.2 

-0.1 
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Figure 27 shows the skill scores averaged over all categories (i.e., AN, 
NN, and BN) by lead time for the Method 1 hindcasts using the NZ SST predictor 
(Table 4). These results show that the hindcast performance tends to decrease 
with increasing lead time, but that there is skill over climatology (LTM based 
forecasts) even at a six month lead time (e.g., HSS > 0 at all leads). 



Oct-Nov Sep-Oct Aug-Sep Jul-Aug Jun-Jul May-Jun 

Lead Times 

- HSS —% Correct - POD - FAR 


Figure 27. Average verification metrics by forecast lead time for hindcasts of 
FlOA PR in Oct-Nov 1970-2009 using Method 1 (tercile matching) 
and the NZ SST predictor. Averaging is done over all forecast 
categories (AN, NN, and BN PR). Verification metrics consist of 
percent correct (% Corr), false alarm rate (FA Rate), probability of 
detection (POD), and Heidke Skill Score (HSS). Refer to Chapter 
II, Section B.6 for more information on the metrics listed in this 
table. Note the general decline in skill as lead time increases, but 
that HSS remains positive at all lead times. 
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2 . 


Composite Analysis Forecasting Method with NZ SST 
Predictor 


Figure 28 shows the hindcast results from applying the composite analysis 
forecast method and the NZ SST to hindcast HOA PR in Oct-Nov at a lead time 
of six months (forecast based on NZ SST in Apr-May). Figure 28a shows the 
historical frequency distributions of PR with respect to NZ SST. The statistically 
significant results are outlined in black and are; all of the BN PR cases, plus the 
AN PR and BN SST cases). The existence of these statistically significant 
results validates the use of the CAF method and indicates that this method 
should perform better than climatology (i.e., better than forecasts based on 
LTMs) (MetEd 2010). 


67 





Figure 28. Results from composite analysis forecasts of Oct-Nov HOA PR 
using the NZ SST predictor at a lead time of six months (i.e., using 
NZ SST in Apr-May). (a) Historical frequency distribution of HOA 
PR probabilities given any SST condition at six months lead. 
Statistically significant results are outlined in black. (b) 
Corresponding probabilistic long-range hindcast of HOA PR in Oct- 
Nov based on AN NZ SST in Apr-May. (c) Corresponding 
probabilistic long-range hindcast of HOA PR in Oct-Nov based on 
BN NZ SST in Apr-May. 


Figure 28 also shows two examples of probabilistic long-range forecasts 
of HOA PR in Oct-Nov when using the CAF Method and the NZ SST predictor 
(when based on the frequency distributions shown in Figure 28a). One forecast 
is based on the occurrence in Apr-May of AN NZ SST (Figure 28b), and the 
other is based on the occurrence in Apr-May of BN NZ SST (Figure 28c). The 
CAF prediction of Oct-Nov HOA PR, given AN NZ SST conditions in Apr-May, 
is: 17% chance of AN PR, 33% chance of NN PR, and 50% chance of BN PR. 
This example forecast indicates a higher (lower) than normal chance of low 
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(high) precipitation. In interpreting these types of forecast results, it is important 
to recall that the normal probability for each of the forecast categories (AN, BN, 
NN) is 33%. So a forecast of a 17% probability indicates a much lower than 
normal probability, and a forecast of a 50% probability indicates a much higher 
than normal probability. A CAF forecast in which one category has a probability 
of 50% or greater indicates that category is more likely than either of the other 
two categories. 

Figure 28a shows a statistically significant result for BN PR given NN SST. 
One would not expect this result, since the predictand-predictor correlation is 
negative, indicating that NN SST would be expected to lead to NN PR. This 
unexpected result may be due to variations in the predictor category as NZ SSTs 
evolve through the northern spring and into the northern summer. During this 
period, the climate system tends to make large adjustments (e.g., onset of Asian 
summer monsoon, initiation of EN and LN events). NZ SST is closely related to 
the SPCZ and to ENLN variations (see Figures 18-24, 26), both of which tend to 
undergo large variations during the northern spring. Thus, unexpected results in 
Figure 28a may represent the effects on the NZ SST predictor of these 
variations. This topic is also discussed in later sections of this report that 
present: (a) results from applying the optimal climate normals approach; and (b) 
hindcast skill when forecasting through the northern spring predictability barrier 
(cf. van den Dool 2007). 

These unexpected results do not show up in the CAF results for shorter 
leads. This indicates that forecasters and forecast end users need to pay extra 
attention to forecast updates to forecasts based on northern spring predictor 
conditions (e.g., updates issued in Jul to forecasts issued in Apr). For the case 
shown in Figure 28, our case, updates to the NZ SST predictor and resulting 
forecasts should be monitored closely as lead time decreases, especially when 
NN SST conditions were observed in the northern spring. 
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3. 


Verification of Hindcasts Based on Aii Years Approach 


Table 5 summarizes the performance of the hindcasts for each of our four 
forecast methods when using NZ SST as the predictor. In this table, the 
verification metrics are averages over all forecast categories (AN< BN< NN) and 
all lead times. All the methods showed skill (e.g., HSS > 0), but Methods 2 and 3 
(linear regression methods) had the highest performance, while Method 1 (tercile 
matching) had the lowest performance. Based on all the verification metrics 
results, we concluded that NZ SST is a viable predictor for generating 
deterministic and probabilistic LRFs. 


Table 5. Average verification metrics for hindcasts of HOA PR in Oct-Nov 
1970-2009 using each of the four forecast methods and the NZ SST 
Predictor. Averaging done over all forecast categories (AN, BN, NN) and 
all lead times of zero to six months. Note the HSS values are all positive, 
indicating forecast skill. See Chapter II, Section B.6 for more information 
on the metrics listed in this table. 



%CoDect 

RJL Rate 

POD 

HSS 

^Method 1 

0.6 

0.5 

0.5 

0.2 

NZ Method 2 

0.8 

0.3 

0.7 

0.5 

laZM^odS 

0.8 

0.3 

0.7 

0.5 

^Method 4 

0.7 

0.5 

0.5 

0.3 


C. LONG-RANGE HINDCASTS: OPTIMAL CLIMATE NORMALS 
APPROACH 

When we examined the time series of the HOA PR predictand and the NZ 
SST predictor, we found the negative correlation between the two variables was 
largely due to the relationship between them during the early and later part of the 
study period, 1970-2009 (Figure 25). When we produced the hindcasts for that 
period using the tercile matching method and the NZ SST predictor, we found 
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that the number of hits was higher at the beginning and the end of the period, 
especially for the last 15 years, 1995-2009. These variable correlation and 
variable skill results indicated that the relationship between the predictor and 
predictand changed during the study period. This may have occurred due to 
changes in the climate system factors that cause variations of the predictor 
and/or predictand. Whatever, the causes, these changes in the predictor- 
predictand relationship indicated that we needed to consider optimal climate 
normals (OCN) approaches to our forecasting (see Chapter II, Section B.5.e). 
OCN approaches allowed us to explicitly account for the changing relationship by 
weighting our forecasts toward the relationships that were dominant in the most 
recent years. The OCN approach is based on the recognition that, in a rapidly 
changing environment, the predictors and forecast methods previously used may 
eventually become non-viable. We applied OCN approaches in the basic and 
advanced forms described in Chapter II, Section B.5.e. 

1. Basic OCN Approach 

The relationship between the HOA PR predictand the NZ SST predictand 
was relatively constant during the last 10-15 years of the 1970-2009 period, as 
indicated by the time series comparisons and the hit rates. Thus, we decided to 
repeat our forecast development and testing process using just the most recent 
15 and the most recent 10 years of data—that is, using data from 1995-2009 
and from 2000-2009 (cf. van den Dool 2007). 

We first did so by taking a basic OCN approach in which we set the tercile 
categories for the predictor and predictand using data from just the most recent 
10 and the most recent 15 years (s opposed to using all 40 years as we did in the 
all years approach (see Chapter III, Section B). This led to changes in the AN, 
BN, and NN categories—for example, the NN HO PR for Oct-Nov category were 
2.9-3.5 mm/day for 1970-2009 and 2.7-3.1 mm/day for 1995-2009. We then 
rehindcasted the most recent 10-year and 15-year periods using the NZ SST 
predictor and the revised tercile categories. 
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To help assess the results of the basic OCN approach, we calculated the 
correlations between the NZ SST predictor and the HOA PR predictand for the 
most recent 10 and 15 years, with the predictor leading by zero to six months 
(Figure 29). The significance levels for the different correlations conducted in our 
study are described in Chapter II, Section B.2. Figure 29 shows that the 
correlations are negative at leads of zero to three months, zero at leads of four to 
five months, and positive at leads of five to six months, including a significant 
positive correlation at a lead of six months. This change in the sign of the 
correlation with lead time was not found when using all forty years of the study 
period (Figure 24). 

Figure 29 also shows the FISS values for hindcasts of 1995-2009 and 
2000-2009 based on tercile matching, and the 10- and 15-year tercile categories 
(or bins) for the predictor and predictand, at all lead times (denoted by 10-yr bins 
and 15-yr bins in the figure legend). The corresponding HSS values when using 
the 40-year tercile categories are also shown {40-yr bins in the figure legend). 
The HSS values decrease with lead time and are approximately zero or less than 
zero (indicating no skill or negative skill) at leads of four to six months. 
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Figure 29. Time series by forecast lead time of; (1) correlations between HOA 
PR predictand and NZ SST predictor; and (2) HSS values for tercile 
matching method hindcasts at lead times of zero to six months, and 
using tercile categories based on the most recent 10, 15, and 40 
years, (a) 1995-2009 correlations (green line); (b) 2000-2009 
correlations (red line); (c) 1995-2009 hindcast HSS values (orange 
and light blue lines); and (d) 2000-2009 hindcast HSS values 
(purple and turquoise lines). The bin periods represent the period 
used to form the predictor and predictand AN, NN, and BN 
categories. Note that the correlations shift from negative to 
positive, and the HSS values decrease, as lead time increases. 


The correlations shown in Figure 29 indicate that the NZ SST predictor 
may not be the optimal predictor for the recent 10 to 15 years. For the 40-year 
period overall, this predictor performed well. But, changes in the climate system 
within that 40-year period may make this predictor less suitable for forecasting 
recent and upcoming years. The correlations in Figure 29 also indicate that 
predictor selection, and forecasting methods in general, may need to vary 
according to the forecast period. In particular, the correlations suggest that other 
variables may be better predictors than NZ SST for the recent 10 and 15 years, 
and upcoming years. 
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2. Advanced OCN Approach 


The results of the basic OCN approach led us to consider a more 
advanced OCN approach in which we repeated our composite and correlation 
analyses using just the most recent 15 years of data. From these analyses, we 
identified two new predictors of HOA PR in Oct-Nov: (a) SST east of the 
Philippines (03N-10N, 143E-152E) for leads of zero to three months; and (b) 
SST west of Namibia (19S-25S, 03W-02E) for leads of four to six months. For 
each of these predictors, there was a negative correlation with the FlOA PR 
predictand. Figure 30 shows representative correlations between global SST 
and our HOA PR predictand, with a focus on the correlations used to identify the 
new predictors. We termed these the multi-region (MR) SST predictors that we 
used when applying the advanced OCN approach to our four forecasting 
methods. 



Figure 30. Correlations of HOA PR in Oct-Nov with global SST for 1995-2009 
with SST leading by: (left panel) zero months (Oct-Nov) and (right 
panel) four months (Jun-Jul). The yellow outlined boxes represent 
the multi-region SST predictor regions identified as having the 
highest correlations to the HOA PR predictand at specific lead 
times when using the advanced OCN approach. We used the SST 
predictor east of the Philippines (left panel) was used for 
forecasting at leads of zero to three months. We used the SST 
predictor west of Namibia (right panel) for forecasting at leads of 
four to six months. Images created using ESRL (2010). 
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All of our statistical significance results using the MR SST predictors 
showed at least one statistically significant relationship for each lead time and 
predictor condition. Thus, we used the advanced OCN approach with each of 
our four forecasting methods, including the CAF method. 

3. Verification of Hindcasts Based on Advanced OCN Approach 

Table 6 summarizes the results of the hindcasts made using forecast 
Method 1 (tercile matching) based on the MR SST predictors. The overall 
hindcast performance was very good, especially for leads of zero to three 
months, and for the AN categories. The NN results were, overall, an 
improvement on those using a non-OCN approach (Table 4). But there were 
some notable cases in which our performance criteria (Chapter II, Section B.6) 
were not met (i.e., POD greater than .50, POD greater than FAR, and/or HSS of 
0.3 or greater were not met). However, the results from using the advanced 
OCN approach were better when applied to the other forecasting methods. 


75 



Table 6. TVerification results from hindcasts generated using Method 1 
(tercile matching). The hindcasts predicted Oct-Nov HOA PR during 
1995-2009 based on the multi-region SST predictors (advanced OCN 
approach; Figure 39), with the predictor leading by zero to six months. 
Columns A, B, C, and D represent the number of hits, false alarms, 
misses, and correct rejections, respectively, for predictions of AN, BN, and 
NN PR. Verification metrics consist of percent correct (% Corr), false 
alarm rate (FA Rate), probability of detection (POD), and Heidke skill 
score (FISS). Compare these advanced OCN verification metrics results 
to those obtained from non-OCN approaches (Table 4). 



A 

B 

C 

D 

Verification Metrics 


Hits 

FA 

Misses 

Corr. Rej. 

% Corr 

FA Rate 

POD 

HSS 

Results for hindcasts at lead time of zero months 

AN PR 

4 

2 

1 

8 

0.8 

0.3 

0.7 

0.6 

BN PR 

2 

2 

1 

10 

0.8 

0.5 

0.5 

0.4 

NN PR 

3 

2 

4 

6 

0.6 

0.4 

0.6 

0.2 

Results ■ 

or hindcasts at lead time of one month 

AN PR 

4 

0 

1 

10 

0.9 

0.0 

1.0 

0.8 

BN PR 

2 

5 

1 

7 

0.6 

0.7 

0.3 

0.2 

NN PR 

3 

1 

4 

7 

0.7 

0.3 

0.8 

0.3 

Results for hindcasts at lead time of two months 

AN PR 

4 

3 

1 

7 

0.7 

0.4 

0.6 

0.5 

BN PR 

2 

2 

1 

10 

0.8 

0.5 

0.5 

0.4 

NN PR 

3 

1 

4 

7 

0.7 

0.3 

0.8 

0.3 

Results for hindcasts at lead time of three months 

AN PR 

4 

1 

1 

9 

0.9 

0.2 

0.8 

0.7 

BN PR 

2 

2 

1 

10 

0.8 

0.5 

0.5 

0.4 

NN PR 

4 

2 

3 

6 

0.7 

0.3 

0.7 

0.3 

Results forhindcasts at lead time of four months 

AN PR 

4 

1 

1 

9 

0.9 

0.2 

0.8 

0.7 

BN PR 

1 

4 

2 

8 

0.6 

0.8 

0.2 

0.0 

NN PR 

2 

3 

5 

5 

0.5 

0.6 

0.4 

-0.1 

Results for hindcasts at lead time of five months 

AN PR 

4 

1 

1 

9 

0.9 

0.2 

0.8 

0.7 

BN PR 

1 

5 

2 

7 

0.5 

0.8 

0.2 

-0.1 

NN PR 

2 

2 

5 

6 

0.5 

0.5 

0.5 

0.0 

Results forhindcasts at lead time of mont 

IS months 

AN PR 

4 

1 

1 

9 

0.9 

0.2 

0.8 

0.7 

BN PR 

1 

3 

2 

9 

0.7 

0.8 

0.3 

0.1 

NN PR 

5 

1 

2 

7 

0.8 

0.2 

0.8 

0.6 
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However, the results from using the advanced OCN approach were better 
when applied to the other forecasting methods. Table 7 shows the skill scores 
averaged over all categories (AN, BN, NN) and all lead times when the advanced 
OCN approach was applied to each of our four forecasting methods. The 
advanced OCN approach yielded higher hindcast verification metrics for each of 
the four forecast methods than were obtained when using non-OCN approaches 
(Table 5). This indicates that using recent predictors and forecast categories 
based on recent years may be improve forecasts of upcoming years. 


Table 7. Average verification metrics for hindcasts of HOA PR in Oct-Nov 
1995-2009 using each of the four forecast methods and the multi-region 
SST predictors (advanced OCN approach; Figure 30). Averaging done 
over all forecast categories (AN, BN, NN) and all lead times of zero to six 
months. Compare these advanced OCN verification metrics results to 
those obtained from non-OCN approaches (Table 5). 



% Correct F.A. Rate 

POD 

HSS 

M-R Method 1 

0.7 

0.4 

0.6 

0.4 

M-R Method 2 

0.9 

0.2 

0.8 

0.7 

M-R Method 3 

0.9 

0.2 

0.8 

0.7 

M-R Method 4 

0.7 

0.4 

0.6 

0.4 


D. LONG-RANGE HINDCAST CASE STUDIES 

In order to further validate our predictive scheme, we tested each LRF 
method, using the NZ SST and MR SST predictors, for individual years within the 
1970-2009 period. We wanted to see how the LRF methods would perform in 
extreme AN and BN years, and the current year. For this, we chose 1997 and 
2006 (AN), 1998 (BN), and 2009 (current year) to use as our hindcast case 
studies. 
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Each hindcast was simulated as if we were the forecasters at the regional 
forecast office, or military installation, asked to predict the coming October- 
November season’s rainfall in the HOA as input for planning an operation. We 
used as our predictors NZ SST and MR SST. We started by analyzing the Apr- 
May SSTs in each predictor region (e.g., NZ or SW Atlantic for the MR), applied 
the analyzed SSTs to the LRF methods, and compared the forecasted PR to the 
actual PR that occurred in the HOA for the October-November period of the year 
being forecasted. We repeated the process for each lead time, and switched to 
the Philippine SST predictor at leads of zero to three months. 

1. 2009 Test Case 

In 2009, the Apr-May analyzed SSTs in the NZ and SW Atlantic predictor 
regions were 23.875C and 22.958C, respectively (Figure 31). 





22.958 C 


»3 

336 

34S 


Figure 31. Composite mean SST (°C) in Apr-May 2009 averaged over: (top 
right) NZ predictor region and (bottom left) MR predictor region. 
The black (red) box represents the NZ (MR) SST predictor regions 
used for analyzing the SSTs at lead times of six months. The area 
averaged SST for Apr-May 2009 was 23.875 (22.958) °C in the NZ 
(MR) predictor regions. Images created using ESRL (2010). 
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The Apr-May NZ and MR SSTs fell into NN SST bins and thus yielded a 
NN precipitation rate forecast for the HOA in Oct-Nov when applying Method 1 
(tercile matching). For Methods 2 and 3, the analyzed SSTs into the formulated 
regression equations based on data from the appropriate years (1970-2009 for 
the NZ SST, 1995-2009 for the MR SST), and the resulting PR value was binned 
into the appropriate tercile category. The NZ predictor region yielded NN 
precipitation rate values of 3.258 (mm/day) and 3.290 (mm/day), for Method 2 
and 3, respectively. The MR SST predictors yielded lower values than the NZ 
SST predictor, yet within the same tercile category of NN precipitation for Oct- 
Nov. The CAF method produced forecasted probabilities of 42% BN rainfall and 
80% NN rainfall in Oct-Nov for the NZ and MR predictors, respectively, given an 
analyzed NN SST in Apr-May (see Chapter II, Section 6, for a description of how 
probabilistic categorical forecasts from the CAF method were verified). 

The actual Oct-Nov PR in the HOA for 2009 was 2.836 (mm/day) which is 
in the NN tercile for the 40-yr tercile bin ranges, and BN for the 15-year tercile bin 
ranges. Recall that the formulation of the tercile bin ranges is dependent on the 
number of years in the dataset, which is why there is a difference between the 
actual PR terciles even though the PR value is the same. 

Therefore, by using the NZ predictor we would have forecasted a NN PR 
for Oct-Nov 2009 given a NN SST in Apr-May, and thus forecasted correctly. 
Using the MR predictor, we would have forecast a NN PR, yet BN precipitation 
occurred and thus forecasted incorrectly. Figure 32 shows the actual and 
anomaly PR reanalysis fields for Oct-Nov 2009. 
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Figure 32. Composite mean (left panel) and anomaly (right panel) PR 
(mm/day) in Oct-Nov 2009. The black outlined box represents the 
HOA PR predictand region. Note that AN, NN, and BN PR 
anomalies occurred within the predictand region (right panel); 
however the actual area-averaged PR was 2.836 mm/day, or NN. 
Images created using ESRL (2010). 

Note that there were AN, NN, and BN anomalies that occurred within the 
predictand region, and therefore the forecast may have been accurate or 
inaccurate, depending on the area within the HOA predictand region. This 
highlights a limitation of area averaged predictand regions—that they cannot 
explicitly represent sub-regional variations. This limitation can be mitigated if the 
region being averaged is small and/or tends to experience spatially 
homogeneous climate variations. Of course, for LRF users and planning lead 
times, this limitation may not be critical (e.g., users who do need high spatial 
resolution for long lead planning). 

We also hindcasted the Oct-Nov 2009 PR in the HOA at each lead time to 

see how the forecasts compared with varying lead time. Table 8 shows the 

results for the 2009 test case. Accurate forecasts are those for which the 

forecasted and observed categories match (e.g., NN forecasted (gray) matches 

NN observed gray, indicating a hit). The categories are color-coded: pink for AN, 

blue for BN, and gray for NN. So, forecast performance can be quickly 

determined by comparing the observed color with the forecast color. Note that 

the gray NN (blue BN) forecasts from the NZ SST (MR SST) predictor tended to 
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perform well. Forecast consistency form one lead to the next can be determined 
by comparing the forecast colors along a row (e.g., all one color along a row 
indicates high consistency). 

Table 8 shows that the forecast performance was low at the shorter lead 
times for all methods, and at most of the lead times using the CAF methods (note 
the many cells in the table that are colored blue (BN) or pink (AN) instead of gray 
(NN)). Overall, the consistency in the forecasts from one lead to the next was 
greater for the NN forecasts, except for the CAF forecasts at the shortest lead 
times. The NZ SST (non-OCN) predictor forecasted slightly better overall than 
the MR SST (OCN) predictor. 

Note from Table 8 that Methods 1-3 using the NZ SST predictor 
performed well at most lead times, but that all the methods, with both predictors, 
did poorly at the shorter leads when they tended to forecast AN PR. This 
occurred because the 2009 SST predictors switched from mainly NN at the 
longer leads to BN at the shortest leads (not shown), thus leading to AN PR 
forecasts at short leads. In summary, the hindcast results for the 2009 case 
were less than ideal, but better than LTM based forecasts would have been, 
when viewed over all lead times. This was a challenging case for all the 
methods, perhaps in part because of the inhomogeneity of the observed HOA PR 
with the HOA (Figure 32) and the variation of the predictors in the transition from 
longer leads to shorter leads. 
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Table 8. Results from hindcasts of HOA PR (mm/day) for Oct-Nov 2009 for 
each forecast method using each of the four forecast methods, and the 
NZ SST (non-OCN) and MR SST (OCN) predictors, at lead times of zero 
to six months. The forecast methods (1-4) are listed at the left. Lead 
times (top row) decrease from left to right. The SST predictor used (NZ or 
MR) is listed in the second row from the top. The forecasted PR values 
are the numerical quantities shown within the columns for each lead time 
and each predictor. The forecasted tercile categories are shown just 
below each forecasted PR. The actual Oct-Nov HOA PR value and 
categories are located at the bottom left (left (right) column shows 
observed tercile category based on non-OCN (OCN) process). Method 1: 
tercile matching; Method 2: averaged linear regression; Method 3: non- 
averaged linear regression; and Method 4; OAF. 



6 months lead 

5 months lead 

4 months lead 

3 months lead 

2 months lead 

1 month lead 

zero months Iea4 

LRF Method 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

1 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

BELOW 

NORMAL 

BELOW 

ABOVE 

NORMAL 

NORMAL 

NORMAL 

ABOVE 

BELOW 

2 

3.259 

3.045 

3.323 

2.851 

3.497 

2.735 

3.436 

2.704 

3.327 

3.019 

3.416 

3.31 

3.666 

3.516 

Category 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

3 

3.245 

3.062 

3.331 

2.852 

3.525 

2.723 

3.468 

2.717 

3.352 

3.033 

3.443 

3.348 

3.698 

3.589 

Category 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

ABOVE 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

4 

42% 

80% 

42% 

67% 

44% 

75% 

67% 

67% 

52% 

67% 

43% 

67% 

67% 

80% 

Category 

BELOW 

NORMAL 

BELOW 

BELOW 

ABOVE 

NORMAL 

ABOVE 

BELOW 

ABOVE 

NORMAL 

NORMAL 

NORMAL 

ABOVE 

ABOVE 
















ACTUAL PRECIP 

2.836 

2.836 














NORMAL 

BELOW 














2. 1997,1998, and 2006 Test Cases 

In 1997, the HOA received a considerable amount of rainfall during the 
short-rains season, along with devastating impacts to the region’s livelihoods 
(ICRC 2004). The HOA PR was well into the AN category. The flooding of 1997 
was followed by drought conditions in 1998, and devastating floods again in 
2006. We investigated these three cases because of their extreme PRs and 
because their actual PR anomalies were relatively coherent within the predictand 
region (Figure 33), unlike the 2009 case. 
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Figure 33. Composite anomalies of PR (mm/day) during Oct-Nov for three 
different years: (left panel) 1997, (middle panel) 1998, and (right 
panel) 2006. The red outlined boxes represent the HOA PR 
predictand region. The HOA experienced extreme AN or BN PR in 
each of these years: AN in 1997 and 2006, and BN in 1998). 
These years were chosen to test the four forecast methods using 
the non-OCN and OCN-based predictor sets. Images created 
using ESRL (2010). 


The hindcast results are shown in Table 9. This table is the same as Table 8 but 
for three different years. Thus, the instructions for interpreting Table 8 (see 
above) apply to interpreting Table 9. Each LRF method performed well in each 
of the three years, with some exceptions (e.g.. Methods 1 and 4 at the longest 
leads). The MR (OCN) SST based forecasts performed better overall. For 1998, 
each LRF method performed well, with the exception of Method 4 using the MR 
SST predictors. Using the Method 4 forecast in isolation would give decision 
makers inaccurate information. But, when combined with the other LRF 
methods, a forecaster might hedge towards a BN PR forecast, since two out of 
three of the other methods forecasted BN PR. The 2006 LRFs at longer leads 
are inconsistent with each other. But, at leads of three to four months and less, 
the forecasts tend toward AN forecasts, in line with actual conditions. 
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Table 9. Results from hindcasts of HOA PR (mm/day) for Oct-Nov of three 
different years using each of the four forecast methods, and the NZ SST 
(non-OCN) and multi-region (OCN) SST predictors, at lead times of zero 
to six months The years are: (top table) 1997, (middle table) 1998, and 
(bottom table) 2006. The forecast methods (1-4) are listed at the left. 
Lead times (top row) decrease from left to right. The SST predictor used 
(NZ or MR) is listed in the second row from the top. The forecasted PR 
values are the numerical quantities shown within the columns for each 
lead time and each predictor. The forecasted tercile categories are shown 
just below each forecasted PR. The actual Oct-Nov HOA PR value and 
categories are located at the bottom left (left (right) column shows 
observed tercile category based on non-OCN (OCN) process). Method 1: 
tercile matching; Method 2: averaged linear regression; Method 3: non- 
averaged linear regression; and Method 4: OAF. Note that all forecast 
method perform well at lead times of zero to four months, but have lower 
skill at the longer lead times, possibly indicating the effects of the spring 
predictability barrier (Chapter III, Section B.2; van den Dool 2007). The 
OCN approach leads to some improvement in performance, especially at 
the longer leads. 



6 months lead 

5 months lead 

4 months lead 

3 months lead 

2 months lead 

1 month lead 

zero months lead 

LRF Method 

NZ 

M-R 

N2 

M-R 

N2 

M-R 

N2 

M-R 

N2 

M-R 

N2 

M-R 

NZ 

M-R 

1 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

2 

3.265 

4.361 

3.337 

4.481 

3.408 

4.335 

3.5975 

4.457 

3.847 

4.505 

3.906 

4.227 

3.692 

4.27 

Category 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

3 

3.24 

4.092 

3.296 

4.297 

3.378 

4.08 

3.58 

4.338 

3.822 

4.367 

3.878 

4.011 

3.633 

4.106 

Category 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

4 

46% 

80% 

46% 

80% 

46% 

80% 

50% 

80% 

50% 

57% 

67% 

100% 

62% 

67% 

Category 

BELOW 

ABOVE 

BELOW 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 
















ACTUAL PRECIP 

4.714 

4.714 














ABOVE 

ABOVE 















6 months lead 

5 months lead 

4 months lead 

3 months lead 

2 months lead 

1 month lead 

zero months lead 

LRF Method 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

1 

BELOW 

NORMAL 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

2 

3.051 

2.975 

2.844 

2.559 

2.642 

2.272 

2.483 

2.219 

2.547 

2.296 

2.632 

2.329 

2.534 

2.311 

Category 

NORMAL 

NORMAL 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

3 

3.208 

3.021 

2.965 

2.584 

2.658 

2.227 

2.464 

2.172 

2.554 

2.268 

2.637 

2.31 

2.568 

2.288 

Category 

NORMAL 

NORMAL 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

BELOW 

4 

45% 

100% 

42% 

80% 

45% 

100% 

53% 

67% 

100% 

67% 

56% 

67% 

57% 

67% 

Category 

BELOW 

NORMAL 

BN/NN 

NORMAL 

BN/NN 

NORMAL 

BELOW 

NORMAL 

NORMAL 

NORMAL 

BELOW 

NORMAL 

BELOW 

NORMAL 
















ACTUAL PRECIP 

2.404 

2.404 














BELOW 

BELOW 
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6 months lead 

5 months lead 

4 months lead 

3 months lead 

2 months lead 

1 month lead 

zero months lead 

LRF Method 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

NZ 

M-R 

1 

BELOW 

ABOVE 

BELOW 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

2 

3.259 

3.568 

3.015 

3.689 

3.167 

3.573 

3.427 

4.074 

3.447 

3.568 

3.553 

3.584 

3.621 

4.211 

Category 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

3 

2.731 

3.419 

2.93 

3.526 

3.116 

3.428 

3.407 

3.903 

3.428 

3.435 

3.532 

3.445 

3.568 

4.025 

Category 

BELOW 

NORMAL 

NORMAL 

ABOVE 

NORMAL 

NORMAL 

NORMAL 

ABOVE 

NORMAL 

NORMAL 

ABOVE 

NORMAL 

ABOVE 

ABOVE 

4 

55% 

67% 

50% 

75% 

38% 

75% 

56% 

75% 

48% 

50% 

64% 

100% 

58% 

60% 

Category 

BELOW 

NORMAL 

BELOW 

ABOVE 

ABV/BN 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 

ABOVE 
















ACTUAL PRECIP 

4.714 

4.714 














ABOVE 

ABOVE 














These three case studies, in which extreme and spatially homogeneous 
PR occurred in the HOA, provide clear evidence that all four forecast methods, 
and the non-OCN and OCN SST predictors, are capable of producing 
consistently accurate LRFs at leads of one to six months. However, the 
performance is noticeable less at leads of five to six months. The LRFs with 
these leads are issued during the northern spring, suggesting that the poorer 
performance the longer lead times may represent the spring predictability barrier 
(Chapter III, Section B.2; van den Dool 2007). In the NH spring, many climate 
processes (e.g., Asian Monsoon, ENSO, lOD, etc) are resetting. As a result, 
predictions that extend through the northern spring tend to have lower skill. 
Thus, forecasts from the longer lead times should be used with extra caution. 

These case study results suggest several conclusions for forecasters 
based on our forecasting methods that may be applicable to other long-range 
forecasting methods. First, the forecasts from different forecasting methods, and 
from non-OCN and OCN processes, might be best used together, rather than 
separately. Thus, for example, a consensus of all the forecasts might in general 
be better than any individual forecast. Second, forecast trends should be 
monitored as lead time decreases, since important changes in the forecasts can 
occur as leads decrease, and especially as the spring predictability barrier is 
passed. Third, the proximity of the forecasted PR values to the category 
thresholds should be monitored and accounted for in issuing and verifying 

forecasts. For example, account for close-call forecasts, such as a forecasted 
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PR value of 3.47 mm/day, which is a borderline NN PR forecast, with the AN PR 
being 3.50 mm/day or above. Fourth, forecasters and forecast users should be 
familiar with the characteristic climate variations patterns for their region of 
interest (e.g., familiar with the results of composite analyses and historical 
frequency distributions associated with ENLN, lOD, etc.). 

E. HINDCAST VERIFICATION SUMMARY 

Figure 34 shows a summary of the hindcast results for each LRF method 
and each predictor (NZ SST / non-OCN) and MR SST / OCN), according to the 
dataset and predictor used. The skill scores are what we termed total verification 
metrics, meaning skill scores averaged over all forecast categories (AN, BN, NN) 
and all lead times. The total verification metrics for hindcasts that used the NZ 
SST (MR SST) predictor are shown in the four leftmost (rightmost) columns. 
Each column is numbered according to the LRF method (see method key in 
lower right). Each skill metric is represented by a different colored symbol (see 
key on left). A reference key at the bottom right shows the type of method for 
each column. For example: the second column from the left (column 2) shows 
the skill metrics for LRF Method 2 (averaged linear regression applied) using the 
NZ SST predictor; and the second column from the right (column 7) shows the 
skill metrics for LRF Method 3 (non-averaged linear regression( using the MR 
SST predictors. To compare the skill of the hindcasts based on the NZ SST 
predictor to that of the LRFs based on the MR SST predictors, compare the left 
and right columns with the same numbers (e.g., left (right) column 3 shows 
Method 3 metrics for NZ SST (MR SST)). 

Figure 34 shows that all of the MSS values are positive, indicating that 
each predictor-predictand pair has skill. In general, the MR SST (OCN) predictor 
has higher skill than the NZ SST predictor, although LRFs based on each of the 
predictors have higher skill than forecasts using the LTM. 

Method 2 and 3 (linear regression LRF models) have the highest total skill 
scores, with the scores from the MR SST based Methods 2 and 3 being higher 
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than the scores from the NZ SST based Methods 2 and 3. The linear regression 
methods may have higher scores than the tercile matching and CAF methods 
because the linear regression methods predict a discrete value from a relatively 
complex model that can accommodate and exploit relatively subtle differences in 
predictor and predictand values. The tercile matching and CAF methods, on the 
other hand, are effectively simpler and cruder models that can only resolve 
predictor and predictand values into one of three broad categories. Thus, the 
linear regression methods tend to be better at describing the details of predictor- 
predictand relationships than the tercile matching and CAF methods. Of the two 
linear regression methods, the non-averaged performed slightly better. This is 
consistent with the idea that better resolution of the details in predictor-predictand 
relationships leads to better forecast performance. 
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Figure 34. Total verification metrics for hindcasts of HOA PR in Oct-Nov 
averaged over all forecast categories (AN, BN, NN) and all lead 
times of zero to six months. Columns 1-4 (5-8) represent 
verification metrics using the NZ (Multi-Region) SST predictor for 
1970-2009 (1995-2009) with each of the forecast methods listed at 
the bottom right of the figure (e.g., tercile matching, linear 
regression and CAF). For example, column 5 represents the 
verification metrics for Method 1 (tercile matching) using the Multi- 
Region SST predictor. Verification metrics consist of percent 
correct (% Corr), false alarm rate (FA Rate), probability of detection 
(POD), and Heidke Skill Score (MSS). Refer to Chapter II, Section 
B.6 for more information on these metrics. The bold black line 
highlights the FISSs. 


Figures 35-38 compare the time series of hindcasted and observed HOA 
PR for Oct-Nov of 1970-2009 for different forecast methods, predictors, and 
lead times. We used these comparisons to help identify the conditions under 
which the hindcasts did performed well and poorly. For example. Figures 35-36 
clearly show that Method 2 and 3 hindcasts based on the NZ SST predictor: (a) 
had better performance before 1984 and after 1993 than in between; (b) had a 
low bias, especially at long leads; and (c) were better at short leads than long 
leads. 
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Figure 35. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 2 (averaged linear regression) and the NZ 
SST predictor for Oct-Nov 1970-2009 at lead times of zero to six 
months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the periods of relatively 
good model performance before 1984 and after 1993 (see Chapter 
III, Section A.3; Figure 24); and (b) the lower variability in the 
hindcasted PR than the observed PR, especially at longer leads. 
All PR values in mm/day. 
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Figure 36. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 3 (non-averaged linear regression) and the 
NZ SST predictor for Oct-Nov 1970-2009 at lead times of zero to 
six months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the periods of relatively 
good model performance before 1984 and after 1993 (see Chapter 
III, Section A.3; Figure 24); and (b) the lower variability in the 
hindcasted PR than the observed PR, especially at longer leads. 
All PR values in mm/day. 


Figures 37-38 show overall very good model performance using Methods 
2 and 3 and the MR SST predictors, especially at shorter leads. These hindcasts 
also show relatively realistic variability, especially at shorter leads. The MR SST 
predictors do better than the NZ SST predictor in hindcasting the amplitudes of 
the interannual variations than relative to the NZ predictor. However, the MR 
SST based hindcasts have relatively low skill between 1999-2002 and then 
again in 2007. The cause for this is not certain, but 1999-2002 was a period of 
NN PR for the HOA, and the forecast models have lower skill forecasting NN PR 
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events. Methods 2 and 3 with the MR predictor were especially skillful in 
hindcasting the extreme events during 1995-2009, especially the AN PR events 
(e.g., 1997, 2004, and 2006) and at shorter leads (Figures 37-38). 
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Figure 37. Time series of actual FlOA PR (black line) and hindcasted PR 
values using Method 2 (averaged linear regression) and the MR 
SST predictors for Oct-Nov 1995-2009 at lead times of zero to six 
months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the overall very good model 
performance, especially at shorter leads; and (b) the relatively 
realistic variability, especially at shorter leads. All PR values in 
mm/day. 
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Figure 38. Time series of actual HOA PR (black line) and hindcasted PR 
values using Method 3 (non-averaged linear regression) and the 
MR SST predictors for Oct-Nov 1995-2009 at lead times of zero to 
six months. The hindcast PR values are dark red (green) at the 
longest (shortest) lead times. Note: (a) the overall very good model 
performance, especially at shorter leads; and (b) the relatively 
realistic variability, especially at shorter leads variability in the 
hindcasted PR than the observed PR, especially at longer leads. 
All PR values in mm/day. 


Figure 39 shows the correlations between the actual Oct-Nov HOA PR 
and the hindcasted HOA PR values shown in Figures 35-38 by lead time. Note 
that the highest correlations are for Methods 2 and 3 based on the MR SST 
predictors, with the Method 2 hindcasts being slightly better correlated than the 
Method 3 hindcasts. The differences in the correlations between the hindcasts 
involving the MR SST predictors and those involving the NZ SST predictors is 
likely due: (a) the smaller number of years represented the hindcasts involving 
the MR SST predictors; and (b) the exclusion from the hindcasts involving the 
MR SST predictors of the problematic years in the middle of the study period 
(1984-1993; Figures 24, 37-38). The correlations involving the NZ SST 
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predictor are more regular across lead times, perhaps due to the larger number 
of hindcasts involved in these correlations. The largest correlations occur at 
leads of two to three months, not at zero months. This may be due to lags 
between when the predictor starts to change and when the predictand responds 
to that change (e.g., maximum variation in predictor precedes by two to three 
months the maximum variation in the predictand). This is plausible, given that 
the predictors and predictand are widely separated, are in two different 
components of the climate system (ocean and atmosphere), and are very likely 
linked via a complex set of low frequency dynamical processes (e.g., ocean 
temperature advection, tropical Rossby and Kelvin waves associated with ENLN 
and lOD mechanisms). 



NZLR Method 2 -A- NZLR Methods —M-R LR Method 2 — I^M-R LR Methoc 3 


Figure 39. Correlations between the actual Oct-Nov HOA PR and the 
hindcasted HOA PR values shown in Figures 35-38 by lead time. 
Forecast method and predictor indicated by colored lines and 
symbols (see key). 
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F. FORECAST FOR OCTOBER-NOVEMBER 2010 

We produced on 01 June 2010 eight LRFs of HOA PR in October- 
November 2010 (Table 10). These eight LRFs were based on all four LRF 
methods, and on the April-May 2010 reanalyzed values for our two SST 
predictor sets (NZ SST and the MR SST). Thus, our forecasts were six month 
lead LRFs. The predictor SSTs were all NN in Apr-May 2010, and all but one of 
the LRFs predicted NN PR in Oct-Nov 2010. This is consistent with the 
transitioning from EN to LN that was occurring in the tropics during Apr-May 
2010. As of the completion of this report, more recent observed predictor values 
were not available. Since these LRFs were issued during the northern spring 
transition, these LRFs extend through the spring transition barrier. Thus, they 
should be used with extra caution and updated LRFs should be monitored by 
forecast users. 


Table 10. Eight long-range forecasts of FlOA PR categories for Oct-Nov 2010 
that we issued on 01 June 2010 at a six month lead time. The LRF 
methods and predictors for each LRF are shown in the table. The 
predictor values were the reanalyzed Apr-May SSTs from the NZ and MR 
predictor regions. Note that all the LRFs predict NN PR in the HOA for 
Oct-Nov 2010, except for one (Method 4, OAF, using the NZ SST 
predictor). Hindcasts generated using these LRF methods and SST 
predictors have their lowest skill at the longest lead (six months). Thus, 
these LRFs should be used with caution and updated LRFs at shorter 
leads should be closely monitored. 


Predictor 

LRF Method 

Forecast PR at 6 months lead 

NZ Predictor 

1 

Normal 

SST: 23.76C 

2 

Normal 


3 

Normal 


4 

Below Normal 

M-R Predictor 

1 

Normal 

SST: 23.03C 

2 

Normal 


3 

Normal 


4 

Normal 
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G. CLIMATE DYNAMICS BASIS FOR LONG-RANGE FORECASTS 

We conducted initial analyses of the climate dynamics that underlie the 
statistical relationships between our HOA precipitation rate predictand and our 
SST predictors north of New Zealand, east of the Philippines, and west of 
Namibia. We conducted this investigation to: (a) identify potential additional 
predictors and predictands; (b) increase our confidence in applying these 
predictors; (c) improve our understanding of the relevant dynamics; and (d) 
determine what processes need to be well modeled in order to improve the 
statistical and dynamical simulation and prediction of HOA precipitation. To 
develop these dynamical analyses, we constructed composites based on the five 
most extreme AN and BN events for our three SST predictors, and for EN and 
LN events. These events are shown in Tables 11 and 12. 


Table 11. The years in which the five most extreme AN and BN conditions 
occurred during Oct-Nov 1970-2009 for each of our three predictors. 



Predictors 

NZ SST 

Philippine SST 

Namibia SST 

Categories 

AN 

BN 

AN 

BN 

AN 

BN 

Years 

1973 

1977 

1995 

1997 

1996 

1997 

1974 

1978 

1998 

2002 

2001 

2000 

1975 

1982 

2000 

2004 

2002 

2004 

1988 

1991 

2001 

2006 

2003 

2006 

1998 

1994 

2007 

2009 

2005 

2007 
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Table 12. The years in which the five most extreme EN and LN conditions 
occurred during Oct-Nov 1970-2009, as indicated by the MEI. 



EN 

LN 


1972 

1971 


1982 

1973 

Years 

1987 

1974 


1994 

1975 


1997 

1988 


1. NZ SST and HOA Precipitation Teieconnections 

Figure 40 shows the composite SSTAs for the five most extreme BN NZ 
SSTs events during Oct-Nov 1970-2009. Recall that BN NZ SST is associated 
with AN HOA PR during Oct-Nov. Note, in this figure: (a) the large area of 
negative SSTAs in the southwestern Pacific centered on our NZ SST predictor 
region (red box) and underlying the SPCZ; (b) the EN pattern of SSTAs in the 
tropical Pacific; (c) the lOD pattern of SSTAs in the 10; (d) the negative SSTAs in 
the tropical Atlantic; and (e) the strong resemblance to the composite SSTAs 
during AN HOA PR (Figure 19). The corresponding SSTAs for the five most 
extreme AN NZ SST events shows nearly opposite patterns. These results 
indicate that our NZ SST predictor is closely related not just to HOA PR but also 
to SPCZ variations, ENLN, lOD, and SST variations in the Atlantic. These 
indications provide information about how NZ SSTs are dynamically related to 
HOA PR. 
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Figure 40. Conditional composite anomalies of sea surface temperature (°C) 
for the five most extreme BN NZ SST events during Oct-Nov 
1970-2009. Recall that BN NZ SST is associated with EN and 
positive lOD events (Figure 17), and with AN HOA PR during Oct- 
Nov. Red box indicates NZ SST region. Red oval indicates area of 
pronounced negative SSTAs in the tropical Atlantic. Image created 
using ESRL (2010). 

Figure 41 shows the composite SSTAs for extreme EN and LN events 
during Oct-Nov 1970-2009. Note: (a) the strong resemblance to the composite 
SSTAs during AN and BN FlOA PR (Figure 19); and (b) the strong resemblance 
to the composite SSTAs during BN NZ SST events (Figure 19); and (c) the 
negative (positive) SSTAs in the tropical Atlantic during EN (LN). These results 
indicate that our NZ SST predictor is closely related to ENLN variations. 
However, the SSTAs in the tropical Atlantic may be better represented by the NZ 
SST predictor than by an ENLN based predictor (e.g., better than by MEI or 
Nino3.4). 
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Figure 41. Conditional composite anomalies of sea surface temperature (°C) 
for the five most extreme EN events (left panel) and LN events 
(right panel) during Oct-Nov 1970-2009. EN and LN events 
identified by use of MEI. Image created using ESRL (2010). 


Rassmusson and Carpenter (1982), as well as many other ENLN studies, 
have shown a close relationship between ENLN and shifts in the locations of the 
SPCZ and ITCZ. In particular, these studies have shown that the SPCZ tends to 
shift to the east (west) during EN (LN) events, while the ITCZ tends to shift to the 
south (north) during EN (LN) events. Figure 42 shows the circulation anomalies 
associated with AN NZ SST, and their overall similarity to the circulation 
anomalies during LN events. Note in particular the indications of enhanced low 
level convergence into the SPCZ, consistent with the indications of enhanced low 
level convergence, convection, and precipitation in the SPCZ during AN NZ SST 
and BN FlOA PR events (Figures 21, 22, and 18, respectively). Recall too that 
AN NZ SST is associated with LN and negative lOD events (Figure 17), and with 
BN HOA PR during Oct-Nov. These results confirm that the NZ SST predictor is 
closely linked to global scale changes in the tropical ocean and atmosphere, 
extending across the Atlantic, India, and Pacific sectors. 
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Figure 42. (Left panel) 3-D schematic of the anomalous Hadley-Walker 
Circulation (HWC) over the Pacific Ocean corresponding to LN 
events and enhanced convection over the MC. (Right panel) 
Conditional composite anomalies of 850 hPa vector winds (m/s) for 
the five most extreme AN NZ SST events during Oct-Nov 1970- 
2009. Recall that AN NZ SST is associated with LN and negative 
lOD events (Figure 17), and with BN FlOA PR during Oct-Nov. 
Bold black line indicates SPCZ. H (L) indicates location of 
anomalous high (low) low-level pressure system. Red arrows 
schematically show the anomalous low-level flow. Image on left 
adapted from Murphree (2009a). Image on right created using 
ESRL (2010). 


Figures 43 and 44 show the composite anomalies of 200 hPa geopotential 
height associated with AN NZ SST and BN NZ SST, respectively. Both figures 
show pronounced anomalous wave trains extending from the maritime continent 
(MC) sector into the extratropics and then back into the tropics over Africa. The 
southern hemisphere anomalous wave train is particularly clear, especially in the 
BN SST case (the case associated with AN HOA PR; Figure 43). These wave 
trains are very likely initiated by convection anomalies in the tropics, especially in 
the eastern 10 - MC - western Pacific region, as indicated by; (1) their apparent 
emanation from that region; (b) their close association with the corresponding 
convection anomalies (cf. Figures 18, 22); (c) their close association with EN and 
LN, for which similar wave trains have been identified (e.g.. Ford 2000) and for 
which tropical upper tropospheric height anomalies similar to those in Figures 
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43-44 are strongly linked (e.g., height anomalies associated with tropical 
Rossby-Kelvin waves induced by anomalous tropical convection during ENLN 
[e.g., Ford 2000; Vorhees 2006]); and (d) their resemblance to wave trains 
initiated by other tropical convection anomalies (Rosencrans 2006). 
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Figure 43. Conditional composite anomalies of 200 hPa GPFI (gpm) for the 
five most extreme BN NZ SST events during Oct-Nov 1970-2009. 
Recall that BN NZ SST is associated with EN and positive lOD 
events (Figure 17), and with AN FlOA PR during Oct-Nov. The red 
arrows schematically show anomalous wave trains that emanate 
from the maritime continent sector, extend into the extratropcial 
northern and southern hemispheres, and then into the tropics near 
Africa. The arrow directions indicate schematically the direction of 
energy propagation in the anomalous wave trains. The Wet \abe\ in 
the upper left indicates that these anomalies are associated with 
above normal precipitation in the HOA. Image created using ESRL 
( 2010 ). 
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Figure 44. Conditional composite anomalies of 200 hPa GPH (gpm) for the 
five most extreme AN NZ SST events during Oct-Nov 1970-2009. 
Recall that AN NZ SST is associated with LN and negative lOD 
events (Figure 17), and with BN FlOA PR during Oct-Nov. The red 
arrows schematically show anomalous wave trains that emanate 
from the maritime continent sector, extend into the extratropical 
northern and southern hemispheres, and then into the tropics near 
Africa. The arrow directions indicate schematically the direction of 
energy propagation in the anomalous wave trains. The Dry label in 
the upper left indicates that these anomalies are associated with 
below normal precipitation in the HOA. Image created using ESRL 
( 2010 ). 


We also investigated the corresponding 200 hPa geopotential height 
anomalies associated with extreme EN and LN events (not shown). These 
composites also provide evidence of anomalous wave trains similar in sign and 
overall pattern to those in Figures 43-44. However, the EN and LN anomalous 
wave trains tend to: (a) be more zonally oriented in the northern hemisphere; and 
(b) have energy propagation paths that extend back into the tropics to the east of 
Africa. 

Taken together, these dynamical analysis results indicate that the NZ SST 
predictor is part of a set of global climate variations that: (a) involve major 
variations in tropical ocean and atmosphere(e.g., ENLN, lOD); (b) include 
extensions into and from the extratropics; and (c) link Africa and nearby regions 
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of the Atlantic and Indian basins to the rest of the tropics and at least parts of the 
extratropics. They also suggest that the NZ SST predictor may be a better 
predictor of HOA PR than ENLN or lOD indices because the NZ SST predictor is 
a better indicator of: (a) anomalies in the strength and location of the SPCZ; and 
(b) the influence of these SPCZ anomalies on anomalous wave trains that extend 
through the southern hemisphere and into Africa. 

Figure 45 shows the composite anomalies of 850 hPa geopotential height 
associated with BN NZ SST in the Atlantic- African region. Comparison of the 
height anomalies in this figure with those in Figure 31 indicate that the Atlantic 
height anomalies in Figure 45 are part of the anomalous wave trains in the 
shown in Figure 43, and that the anomalous wave train has equivalent barotropic 
structure. The associated wind anomalies (red arrows in Figure 45) indicate that 
the wave trains contribute to anomalous convergence into the HOA, consistent 
with the association of BN NZ SST with AN HOA PR. 
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Figure 45. Conditional composite anomalies of 850 hPa GPH (gpm) for the 
five most extreme BN NZ SST events during Oct-Nov 1970-2009. 
Recall that BN NZ SST is associated with EN and positive lOD 
events (Figure 17), and with AN FlOA PR during Oct-Nov. H (L) 
indicates location of anomalous high (low) low level geopotential 
heights. The red arrows show schematically the anomalous wind 
directions. Note the anomalous convergence into the HOA region. 
Image created using ESRL (2010). 


Figure 46 shows the composite anomalies of precipitation rate (PR) 
associated with BN NZ SST in the Atlantic-African region. Note: (a) the positive 
PR anomalies over the southern Ethiopia, Somalia, and west-central Kenya 
portion of the HOA, consistent with the corresponding anomalous convergence 
into that region shown in Figure 45; and (b) the negative PR anomalies to the 
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east of the HOA, consistent with the corresponding anomalous high pressure 
over that region shown in Figure 45 and the negative tropical Atlantic SSTs 
shown in Figure 40. 



Figure 46. Conditional composite anomalies of precipitation rate (PR; mm/day) 
for the five most extreme BN NZ SST events during Oct-Nov 
1970-2009. Recall that BN NZ SST is associated with EN and 
positive lOD events (Figure 17), and with AN FlOA PR during Oct- 
Nov. The red box indicates our HOA PR predictand region. Image 
created using ESRL (2010). 


The results in Figures 45-46 indicate that wave trains associated with BN 
NZ SST contribute to circulation variations over Africa that lead to increased 
precipitation in the HOA. Similar analyses of conditional composites for extreme 
AN NZ SST events show approximately opposite patterns. 


104 












2 . 


MR SST and HOA Precipitation Teieconnections 


The teieconnections between the MR predictors and the HOA are similar 
to the teieconnections between ENSO and lOD and the HOA, although there are 
some amplified features that are worth noting. 

a. Teieconnections at 4-6 Month Lead Time 

At leads of four to six months, we found that SST west of Namibia 
was a useful predictor of HOA PT in Oct-Nov (Chapter III, Section C.2). In this 
region a key feature of the climate system is the South Atlantic Subtropical High 
(SASH). Several prior studies have linked anomalies in the structure and 
intensity of the SASH have been dynamically linked to changes in east African 
rainfall (e.g., Nicholson and Entekhabi 1987; Kabanda and Jury 1999; Mutai and 
Ward 2000; Wang et al. 2009). 

The normal low-level circulation pattern around the SASH causes 
southerly flow along the southwestern coast of Africa and easterly trades along 
the equator, both of which cause upwelling and LTM SSTs that are relatively cool 
in those regions. In the center of the SASH, LTM SSTs are relatively warm due 
to low wind speeds and convergent Ekman transport into that region. Figure 47 
shows that four to six months prior to the AN HOA PR in Oct-Nov, the 850 hPa 
geopotential heights tend to be higher than normal on the southeastern side of 
the LTM position of the SASH, with corresponding southerly and easterly wind 
anomalies on the eastern and northern flanks of the SASH that are stronger than 
normal winds along the coast of southwest Africa and across much of the 
subtropical South Atlantic (red arrow in Figure 47). 
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Figure 47. Conditional composite anomalies of 850 hPa GPH (gpm) for the 
five most extreme BN Namibia SST events during Apr-May 1970- 
2009. Recall that BN Namibia SST is associated with AN HOA PR 
during Oct-Nov. The H indicates an anomalous high-pressure 
region. The direction of the corresponding wind anomalies is 
shown schematically by the red arrow. The cooler SSTs region is 
explained in the main text discussion of Figure 40. Image created 
using ESRL (2010). 


Figure 48 shows the SST anomalies corresponding to the height 
anomalies shown in Figure 47. Note the regions of negative SSTAs occurring 
under and to the northeast of the southeasterly wind anomalies marked in 
Figures 47 and 48. The negative SSTA region coincides with our Namibia SST 
predictor region (Chapter III, Section C.2). These SSTAs are consistent with the 
wind driven heat flux anomalies and upper ocean Ekman transport anomalies 
that would be expected from the southeasterly wind anomalies. This suggests 
that the Namibia SST predictor is associated with anomalies in the SASFI. The 
results shown in Figures 19-20, 22, and 43-45 suggest that the SASFI 
anomalies may be driven by anomalous wave trains that: (a) are initiated by 
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anomalous SST and convection in the eastern 10—western Pacific region, 
including anomalies in the SPCZ region; and (b) extend through the southern 
hemisphere into the South Atlantic. 



Figure 48. Conditional composite anomalies of SST (°C) for the five most 
extreme BN Namibia SST events during Apr-May 1970-2009. 
Recall that BN Namibia SST is associated with AN HOA PR during 
Oct-Nov. The negative SSTA region west of Namibia corresponds 
to the Namibia SST predictor region identified as having the highest 
correlations, at leads of four to six months, to the HOA PR 
predictand for 1995-2009 (see Chapter III, Section C.2). The H 
indicates an anomalous lower tropospheric high-pressure region 
and the black arrow indicates the corresponding wind anomaly (see 
Figure 47). Image created using ESRL (2010). 


b. Teleconnections at 0-3 Month Lead Time 

At leads of zero to three months, we found that SST east of the 
Philippines was a useful predictor of HOA PT in Oct-Nov (Chapter III, Section 
C.2). Figure 49 shows the correlations between HOA PR in Oct-Nov and global 
SST in Oct-Nov (zero lead) for our full 40-year study period, 1070-2009, and our 
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15-year OCN period. The two correlation maps are very similar but the 
correlations are stronger for the shorter period. In particular, the negative 
correlations east of the Philippines are much stronger for the shorter period and 
are the strongest correlations globally for this short period. 



Figure 49. Correlations of HOA PR in Oct-Nov with global SST for 1970-2009 
(left panel) and 1995-2009 (right panel), with SSTs leading by zero 
months (Oct-Nov). The yellow outlined boxes represent NZ 
predictor region (left panel) and the Philippine SST predictor region 
(right panel). Note the similar patterns for the two periods, but with 
stronger correlations for 1995-2009. Images created using ESRL 
( 2010 ). 

The importance of the Philippine SST predictor, and its greater 
importance in the most recent 15 years, may be related to the long term changes 
in SST that have occurred in this region and elsewhere. Figure 50 shows the 
global SSTAs for Oct-Nov of 1970-2009. The LTM base period for the 
calculation of these anomalies is 1968-1996 (see Chapter II, Section B.2), so 
these anomalies represent decadal scale changes in SST (e.g., changes 
associated with global warming). Note in particular that the SSTs in the 
Philippine SST predictor region and the NZ SST predictor region increased by 
approximately 0.5° C compared to the LTM values. AN SST in these predictor 
regions during Apr-Nov is associated with BN FlOA PR in Oct-Nov. Thus, the 
decadal scale SST increases in these two regions might be expected to lead to 
more BN HOA PR events. The HOA has experienced a higher frequency of 
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recurring droughts and floods in the last decade relative to the last 40 years (e.g., 
HOA PR was below its LTM value in nine out of the 15 years during 1995-2009). 



Figure 50. Composite anomalies of SST (°C) for Oct-Nov of 1995-2009. Note 
the positive SSTAs, especially in the northern hemisphere and in 
the western North and South Pacific. Image created using ESRL 
( 2010 ). 

3. Summary of Climate Dynamics Analyses 

Based on our initial analyses of the climate dynamics, we developed 
schematic representations of the integrated set of global scale patterns and 
processes associated with AN and BN PR in the HOA in Oct-Nov (Figures 51- 
52). These schematic figures describe lower tropospheric and sea surface 
conditions, but the low level atmospheric circulation features are due in part to 
anomalous wave trains that extend throughout the troposphere. These 
schematics primarily describe the conditions associated with AN and BN periods 
for our NZ SST predictor. However, similar schematics apply to AN and BN 
periods for our Philippine and Namibia SST predictors. 

Figure 51 shows the conditions associated with AN HOA PR and BN NZ 
SST. These are the conditions in the BN NZ SST is likely to perform well in 
predicting AN PR in our HOA predictand region. Figure 52 shows the conditions 
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associated with BN HOA PR and AN NZ SST. These are the conditions in which 


the AN NZ SST is likely to perform well in predicting BN PR in our HOA 
predictand region. 
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Figure 51. Schematic diagram of lower tropospheric and sea surface 
conditions that occur in October-November when: (a) precipitation 
is greater than normal in the HOA (i.e., during AN HOA PR events); 
and (b) when SSTs north of New Zealand are cooler than normal 
(i.e., during BN NZ SST events) and. Map background from 
Google Maps (www.maps.google.com). 
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Figure 52. Schematic diagram of lower tropospheric and sea surface 
conditions that occur in October-November when: (a) precipitation 
is less than normal in the HOA (i.e., during BN HOA PR events); 
and (b) when SSTs north of New Zealand are warmer than normal 
(i.e., during AN NZ SST events) and. Map background from 
Google Maps (www.maps.google.com). 


Figure 53 shows the conditions associated with above normal precipitation 
for the greater FlOA region. These conditions are very similar to those in Figure 
51, but there are some notable differences in the circulation and precipitation 
patterns over and near Africa. Note that the AN PR region is located further to 
the north than in Figure 51. The differences between Figures 51 and 53 occur 
because Figure 51 is focused on conditions associated with BN NZ SST events, 
and these events are focused on the southern and central portions of the greater 
HOA region, due to our choice of predictand region (se Chapter II, Section B.1) 
for which we then identify NZ SST as a good predictor. As discussed in Chapter 
II, the selection of the predictand is critical in determining the outcomes from the 
statistical and dynamical analyses, and the selection of predictors, the 

development of forecast methods, and the results from hindcasts and forecasts. 
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Figure 53. Schematic diagram of lower tropospheric and sea surface 
conditions that occur in October-November when precipitation is 
greater than normal in the HOA (i.e., during AN HOA PR events). 
Map background from Google Maps (www.maps.google.com). 

In summary, the major patterns associated with anomalously high (low) 
precipitation in the HOA during Oct-Nov: 

1. Anomalously cool (warm) SSTs in the western tropical Pacific, 
including the regions north of New Zealand and east of the Philippines, 
the eastern tropical Indian, and the tropical Atlantic basins 

2. Anomalously warm (cool) SSTs in the central and eastern tropical 
Pacific and the western tropical Indian basins 

3. Anomalously weak (strong) easterlies over the central and eastern 
tropical Pacific basin 

4. Anomalously strong (weak) easterlies over the western tropical Indian 
basin 
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5. Increased (decreased) sea level pressure (SLP) over the maritime 
continent (MC) and decreased (increased) SLP over eastern Africa 

6. Anomalous low level convergence (divergence) into eastern Africa, 
associated in part with anomalous wave trains emanating from the MC 
sector and extending through the extratropics and back into the tropics 
near Africa 

7. Anomalies in the Walker circulation that lead to increased (decreased) 
upward motion and convection over the HOA and decreased 
(increased) upward motion and convection over the MC 
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IV. SUMMARY, CONCLUSION AND RECOMMENDATIONS 


A. ANALYSIS OF FORECASTING METHODS AND RESULTS 

This study explored the viability of using state-of-the-science climate 
dataset and methods to develop skillful LRFs of October-November rainfall in the 
HOA of leads of up to six months. The primary aim of this research was to 
generate LRF support for military and non-military operations in the HOA. We 
highlighted that extreme climate events may act as a catalyst for creating or 
exacerbating the need for operations in the HOA, many of which rely on accurate 
climate forecasts for their decision making and environmental awareness. We 
chose to focus our research on the HOA because of its interest to military and 
non-military operations alike, and because little has been done in regards to 
developing statistically based precipitation forecasts for the HOA out to six 
months. 

We employed several analysis methods (outlined in Chapter II) using 
readily available state-of-the-science reanalysis and datasets to conduct 
analyses of the LTM and climate variations in HOA rainfall. The analyses 
allowed us to identify global scale interannual patterns and processes related to 
variations in HOA rainfall, including several climate teleconnections that had not 
previously been identified and/or investigated for their long-range forecasting 
potential. 

Based on established criteria outlined in Chapter I, we chose a 
customized predictand region within the HOA for determining the potential for the 
long-range forecasting of precipitation. After correlating our predictand with 
several climate variables and indices, we found that SSTs north of New Zealand 
(NZ) had the highest significant correlations out to six-month lead times. In 
addition, we applied the optimal climate normals (OCN) concept for long range 
forecasting, using two different approaches. From our OCN based correlation 
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analyses, we identified two additional predictors of HOA precipitation in Oct-Nov: 
(1) SSTs east of the Philippines for leads of zero to three months; and (2) SST 
west of Namibia for leads of four to six months. 

We then used the NZ SST predictor to generate hindcasts of HOA 
precipitation for Oct-Nov 1970-2009 at leads of zero to six months four different 
deterministic and probabilistic forecasting methods (e.g., tercile matching, linear 
regression, composite analysis forecasting). We also applied the same 
forecasting methods using the Philippine and Namibia SST predictors to hindcast 
HOA precipitation for Oct-Nov 1995-2009 at leads of zero to six months. We 
also conducted and analyzed LRFs for five specific cases to: (a) better determine 
the performance of our predictors and forecast methods; and (b) issue an 
experimental forecast of precipitation in the HOA for Oct-Nov 2010. 

Our results indicate that skillful LRFs of HOA Oct-Nov precipitation are 
possible out to six month leads for the HOA using our specifically identified SST 
predictor regions in the Pacific and Atlantic. The skill when using these 
predictors exceeds that from using more established climate indices as 
predictors (e.g., ENLN and lOD indices). 

We concluded by summarizing the climate dynamics that underlie the 
predictor-predictand relationships that we used to generate our hindcasts and 
forecasts. 

B. RECOMMENDATIONS 

It is clear that climate analyses and long-range forecasts based on the 
advanced climate datasets and methods used in this study can provide 
immediate support to the improvement of military and non-military plans and 
operations in the HOA. The following section highlights some of the key areas of 
future research that can support these improvements. 

1. Our preliminary results were encouraging, but further verification is 
needed by testing the LRF methods in operational settings and exploring our key 
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findings in the context of long-range planning for DoD and non DoD operations 
(e.g., FEWS NET, AFRICOM, 14 WS, etc). Military planners generally fuse 
climate data products into the planning cycle several months to several years in 
advance of any possible contingency (Hernandez and Wunder 2004). The LRF 
processes and methods developed in this study are unique in that they provide 
planners the ability to adjust their contingency plans for the HOA based on 
specific climate variations (e.g., EN, positive phase of the lOD, etc.). Using 
analyses of the regional impacts of extreme climate events, planners can link the 
effects of each climate variation to the region’s socio-economic and political 
sectors that induce specific military and non-military actions. For example, when 
LRFs predict AN Oct-Nov rainfall in the HOA, and therefore, likely flooding, 
military operations may be oriented more towards humanitarian activities such as 
improving roads and bridges, and hardening shelters for local communities 
several months before the short rains season. If the LRFs predict a dry short 
season, the military may help communities prepare by improving water supply 
systems or by providing peacekeeping and security operations where dry 
conditions are likely to lead to border conflicts or civil unrest. 

2. Expand the type of research we conducted to additional HOA 
seasons, climate variations, and sub-regions in an effort to fit better with existing 
forecast products (e.g., products from Kenyan Meteorological Department, 14 
WS). 

3. Explore inter-decadal and decadal climate variations that may 
influence PR in the HOA. The decadal scale variations in our correlations in the 
skill of our hindcasts indicate that more research is needed on the causes of 
lower frequency variations and their impacts on the HOA. This research will in 
turn support improvements in long-range forecasting for the HOA. 

4. Apply the methods to other operational regions. The readily 
accessible advanced climate datasets and straightforward methods we used in 
this study can be adapted to other regions of the world and various operations. 
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5. Develop a user-friendly online application that produces the 
deterministic and probabilistic forecasts generated in this study. The data 
manipulation needed to conduct our analyses and produce our forecasts was 
relatively complex and lengthy, but it could be automated and made much easier 
and more accurate if it could be done via a well designed and tested online 
application. 

6. Expand the consideration of partnerships between DoD and civilian 
organizations, who do advanced climate research and operational climate 
support (e.g., NOAA), to share ideas on climate research, and climate analyses 
and forecast products. This could be especially useful in helping DoD leverage 
existing resources and expertise to rapidly improve support for DoD planning and 
operations. 
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